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ABSTRACT 


The  PHALANX  Close-In-Weapons  System  provides  the  last  layer  of  defense 
against  an  incoming  anti-ship  missile.  Current  dis|>ersion  data  shows  that  the 
rotating  six-barrel  system  firing  at  75  shots/second  gives  variable  dispersion 
greater  than  2  mradian.  This  results  in  a  4  m  spread  at  its  maximum  2000  m 
range.  Computer  simulations  have  shown  that  a  fast  missile  can  penetrate  to 
within  a  couple  hundred  yards  of  the  ship,  where  even  if  it  was  destroyed,  the 
resulting  ballistic  debris  would  hit  the  ship  in  a  shotgun  blast  pattern. 


In  order  to  investigate  the  physics  behind  the  dispersion  and  test  future 
design  changes,  a  detailed  finite  element  model  of  the  PHALANX  gun  was 
developed.  This  is  the  first  finite  element  model  of  the  gun  ever.  The  normal 
modes  of  a  single  barrel  have  been  matched  favorably  with  those  of  a  real  barrel. 
This  model  has  been  used  to  calculate  the  normal  modes  of  vibration  of  the  whole 
gun  in  order  to  predict  the  motion  of  the  complex  system.  Barrel  tip  displacement, 
which  plays  a  critical  role  in  the  dispersion  pattern,  was  determined  after  applying 
the  periodic  firing  excitation  force  and  damping.  Its  motion  corresponds  well  to  the 
2  mradian  dispersion.  In  addition,  its  first  mode  of  vibration  at  14  Hz  appears  to 
be  significant  in  the  analysis  of  the  dispersion  data.  This  analysis  has  focused 
attention  on  certain  components  of  the  gun,  such  as  the  double  angular  contact 
bearing,  which  may  be  a  leading  cause  of  the  dispersion  pattern.  Finally,  this  gun 
model  can  also  support  design  changes  such  as  the  addition  of  barrel  restraints  or 


the  proposed  new  barrels  and  their  effect  on  dispersion  predicted  prior  to  any 


actual  hardware  tests  or  procurement  orders. 


iii 


Table  of  Contents 


I.  INTRODUCTION  .  1 

II.  MOTIVATION  FOR  STUDYING  THE  PHALANX  GUN .  5 

A.  DESCRIPTION  OF  THE  PHALANX  SYSTEM .  5 

B.  PHALANX  ROUND  DISPERSION .  10 

C.  THE  BALUSTIC  DEBRIS  PROBLEM .  15 

III.  FINITE  ELEMENT  THEORY .  22 

A.  THE  DECISION  TO  USE  FINITE  ELEMENT 

MODELING .  22 

B.  THE  DISCRETIZATION  OF  A  SIMPLE  ROD .  25 

C.  MATRDC  FORM  OF  EQUILIBRIUM  EQUATIONS  .  30 

D.  FREE  VIBRATION  IN  A  STRUCTURE  OF  ELEMENTS .  34 

IV.  MODEL  DESCRIPTION  .  41 

A.  TYPES  OF  FINITE  ELEMENTS  USED  IN  THE  MODEL  .  41 

1.  Brick  Elements .  42 

2.  Thin  Shell  Elements  .  44 

3.  Beam  Elements  .  45 

4.  Spring  Elements .  46 

5.  Rigid  Elements  .  47 

B.  CONSTRUCTION  OF  THE  PHALANX  GUN  .  49 

1.  Gun  Barrels  .  50 

2.  Muzzle  Clamp .  50 

iv 


3.  Mid-Barrel  Clamp 

4.  Stub  Rotor . 


51 


52 

5.  Rotor .  52 

6.  Double  Row  Contact  Bearing .  54 

7.  Needle  Bearing .  57 

8.  Gun  Body  .  60 

9.  Recoil  Adapters  .  61 

10.  Ball  Joint  .  63 

C.  BOUNDARY  CONDITIONS  .  65 

D.  COMPARISON  BETWEEN  THE  REAL  AND  MODELED 

GUN  68 

V.  MODES  OF  VIBRATION .  71 

A.  MODES  OF  VIBRATION  OF  THE  GUN  BARRELS  .  71 

B.  MODES  OF  VIBRATION  OF  THE  COMPLETE  GUN 

MODEL  .  75 

1.  First  Flexible  Mode .  78 

2.  Second  Flexible  Mode  .  80 

3.  Third  Flexible  Mode .  83 

4.  Fourth  Flexible  Mode  .  84 

5.  Fifth  Flexible  Mode .  85 

6.  Sixth  Flexible  Mode  .  86 

C.  MODAL  TRUNCATION .  86 

VI.  DRIVING  AND  DAMPING  OF  THE  PHALANX  GUN  .  88 

A.  DYNAMIC  ANALYSIS  OF  THE  GUN  MODEL .  88 

B.  DAMPING  IN  THE  PHALANX  GUN .  90 


V 


C.  DRIVING  THE  PHALANX  GUN 


96 


D.  RESPONSE  OF  THE  GUN  MODEL .  104 

E.  FREQUENCY  RESPONSE  FUNCTIONS .  107 

F.  OTHER  FACTORS  AFFECTING  DISPERSION  .  Ill 

VII.  SUGGESTIONS  FOR  FUTURE  WORK .  115 

A.  REFINING  THE  GUN  MODEL .  1 15 

1.  Muzzle  and  Mid-Barrel  Clamps  .  1 16 

2.  Gun  Body  Housing .  117 

3.  Rotor .  119 

4.  Stub  Rotor .  120 

5.  Recoil  Adapters  .  120 

B.  VAUDATION  OF  THE  MODEL  THROUGH  TESTING  .  121 

C.  ADDITIONAL  MODELING .  122 

1.  Modifying  the  Existing  Firing  Force .  122 

2.  Adding  Additional  Forces .  123 

3.  Additional  Modeling  .  123 

VIII.  CONCLUSIONS .  126 

APPENDIX  A.  GUN  BARREL  NORMAL  MODES .  129 

APPENDIX  B.  GUN  MODEL  MODE  SHAPES  .  134 

LIST  OF  REFEPJINCES .  173 

INITIAL  DISTRIBUTION  UST .  176 


VI 


ACKNOWLEDGEMENT 


The  author  is  grateful  for  the  xinbounded  guidance  and  support  provided  by 
Dr.  W.  B.  Colson  and  Dr.  J.  H.  Gordis.  The  author  would  also  like  to  thank  Mike 
Hatch  for  his  volunteerism  and  enthusiasm  in  making  this  project  all  that  much 
better.  Also,  the  author  would  like  to  thank  his  partner,  LT  Don  MacNeU,  who 
provided  the  drive  in  helping  finish  this  work  in  such  a  short  time.  In  addition,  the 
author  especially  would  hke  to  thank  his  parents.  Bill  and  Sue,  who  have  given  the 
author  more  than  he  could  ever  have  asked  for.  Thanks  also  goes  to  the  author’s 
lifelong  friends,  Jon  and  Pam  Benson,  Bill  Knowles  and  Joe  Murphy  for  all  their 
encoiiragement.  Finally,  the  author  wishes  to  give  a  special  thanks  to  Pattie 
McCluskey  who  inspired  the  author  during  the  final  mile. 


vii 


I.  INTRODUCTION 


The  PHALANX  Close-In-Weapon-System  (CIWS)  is  considered  to  be  a  ship’s 
last  layer  of  defense  against  incoming  missiles.  The  PHALANX  CIWS  is  shown  in 
Figure  1-1.  Nearly  every  ship  of  the  current  465  ship  U.  S.  fleet  has  at  least  one  of 
these  gun  mounts.  For  some  ships,  it  is  the  only  means  of  missile  defense. 
Although  it  has  been  in  the  fleet  for  over  a  decade,  it  has  never  been  tested  in 
combat,  so  its  efiectiveness  has  largely  been  measured  by  either  missile  range  tests 
or  computer  simiilations.  Within  the  past  flve  years,  there  has  been  documented 
concern  that  this  defensive  weapon  system  may  not  be  that  effective,  especially  in 
defeating  sophisticated  missiles  that  are  currently  available  to  Third  World 
countries.  The  problem  of  its  effectiveness  is  attributed  almost  entirely  to  that  of 
the  gun  itself  and  not  its  Are  control  system  [Ref.  1].  At  the  maximum  range  of  the 
PHALANX  gun,  the  incoming  missile  or  patrd  boat  presents  a  small  cross-section. 
The  angular  dispersion  of  bullets  leaving  the  gun  can  dimmish  the  gun’s  ability  to 
hit  targets,  especially  as  the  target’s  size  decreases,  becomes  more  maneuverable, 
or  increases  its  velocity. 

As  part  of  a  class  project,  supervised  by  Professor  W.  B.  Colson,  the  problems 
associated  with  the  gxm  were  illustrated  in  some  of  the  best  case  scenarios. 
Specifically,  a  computer  simulation  was  used  to  study  a  PHALANX  engagement  of 
an  incoming  anti-ship  missile.  The  results  of  this  class  project,  discussed  in 
Chapter  II,  provided  the  motivation  for  this  thesis. 

The  purpose  of  this  study  has  been  to  use  a  sophisticated  computer  code  to 
create  a  finite  element  model  of  the  PHALANX  gun  in  order  to  better  imderstand 
its  complex  dyneimic  behavior.  The  Naval  Postgraduate  School  has  access  to  one  of 
the  most  advanced  finite  element  codes  available  in  private  industry  today,  called 
I-DEAS  (Integrated  Design  Engineering  Analysis)  which  is  developed  and 
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Figure  1*1.  The  PHALANX  Cloee-ln  Weapon  System  [Ref.  2], 


supported  by  SDRC  (Structural  Dynamics  Research  Corporation).  The  finite 
element  model  was  used  to  predict  barrel  tip  displacement  through  such  causes  as 
barrel  oscillations  and  the  coupling  of  the  barrels  to  the  gun  mount  through  its 
bearings.  In  addition,  the  finite  element  model  provides  an  important  tool  for  the 
evaluation  of  future  design  changes  of  the  PHALANX  gun  system.  Chapter  III 
discusses  the  benefits  of  using  this  type  of  vomputer  code  and  the  theory  behind  its 
development. 

Chapter  IV  then  discusses  how  raw  data  from  blueprints  and  technical 
discussions  with  engineers  familiar  with  the  PHALANX  CIWS  was  molded  into  a 
model  of  the  gxm.  Field  trips  to  Hughes  Missile  Systems  Company,  Pomona, 


California,  and  the  USS  ARKANSAS  provided  further  modeling  information. 
Valuable  finite  element  modeling  experience  by  Mike  Hatch,  an  expert  on  rotating 
computer  disk  systems  >ind  Professor  J.  H.  Gordis,  a  former  SDRC  mechanical 
vibration  theorist,  ir^atly  enhanced  the  effectiveness  of  the  model. 

In  Chapter  V,  the  normal  modes  of  vibration  are  calculated  for  the  PHALANX 
gun.  Solving  for  the  normal  modes  is  an  important  and  insightful  step  into 
understanding  the  dynamic  response  of  any  structure.  Determining  which  normal 
modes  of  a  structure  are  excited  after  applying  particular  excitation  forces  and 
setting  boundary  conditions  ultimately  determines  how  that  structure  will  behave. 
In  analogy,  any  color  of  light  can  be  made  up  of  the  three  primary  colors,  red,  blue, 
and  yellow.  The  response  of  any  system  can  be  defined  in  terms  of  its  normal 
modes.  The  dynamic  response  package  of  I-DEAS  uses  modal  degrees  of  freedom 
instead  of  physical  degrees  of  fineedom,  in  order  to  greatly  reduce  the  dynamic 
response  solving  time. 

Chapter  VI  discusses  how  these  normal  modes  are  used  in  a  dynamic 
response  analysis.  Damping  is  defined  for  the  gun  system  and  the  primary  driving 
force  of  the  gun  while  firing  is  applied.  The  resulting  motion  of  the  barrel  tip, 
which  gives  a  good  indication  of  the  dispersion  pattern,  is  then  displayed  for  each 
of  the  three  translational  directions.  More  importantly,  the  firequency  response 
functions  are  also  calculated  and  are  used  to  reveal  which  normal  modes  are  being 
excited  and  to  what  extent.  Such  an  analysis  may  ultimately  lead  to  tuning  a 
problem  mode  so  that  its  excitation  is  limited,  and  the  resulting  barrel  tip  motion 
is  correspondingly  reduced. 

Recommendations  for  follow-on  modifications  of  the  model  are  listed  in 
Chapter  VII.  Although  the  current  model  has  been  correlated  in  a  number  of  ways 
with  the  real  gun,  it  is  important  to  refine  the  model  further  so  as  to  enhance  its 
accuracy  while  reducing  its  total  number  of  degrees  of  fireedom.  By  reducing  the 
model’s  degrees  of  freedom,  the  computer  solve  time  is  shortened  which  can  lead  to 
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a  smoother  animation  picture,  thus  providing  better  insight  into  its  dynamic 
behavior.  Also,  if  the  proposed  muzzle  restraint  and  isogrid  is  added,  it  will  be 
important  to  keep  the  model’s  complexity  within  a  manageable  budget. 

Although  the  primary  purpose  of  this  study  was  to  huild  and  demonstrate  the 
practical  value  of  a  finite  element  model  of  the  gun,  some  preliminary  analysis  of  it 
was  performed  and  the  results  are  discussed  in  Chapter  VIII.  Of  key  interest  is 
the  first  flexible  mode  of  vibration  at  14  Hz  and  the  double  angular  contact  bearing 
and  the  effects  they  may  have  on  the  dispersion  pattern. 

The  purpose  of  this  study  has  been  multi-faceted.  While  this  computer  model 
will  be  important  in  predicting  the  observed  dispersion  pattern,  it  can  also  serve  as 
a  program  manager’s  tool  in  helping  make  future  design  decisions.  Components 
can  be  swapped  out  with  new  ones,  and  their  effect  on  the  dynamic  behavior  of  the 
gun  can  be  predicted  long  before  any  hardware  tests  are  actually  performed.  In 
summary,  a  tool  of  the  1990’s  is  being  used  to  improve  1950’8  technology. 
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II.  MOTIVATION  FOR  STUDYING  THE  PHALANX  GUN 


A.  DESCRIPTION  OF  THE  PHALANX  SYSTEM 

The  PHALANX  Close-In  Weapon  System  (CIWS)  is  a  ship-based,  stand  alone 
gun  mount  designed  to  disable  an  incoming  missile  or,  in  a  recent  application, 
engaging  fast  patrol  boats.  It  has  been  actively  deployed  for  over  a  decade  in 
virtually  every  ship  of  the  U.  S.  fleet,  and  has  had  to  undergo  continual  review  and 
various  modifications  to  match  the  ever  ’’fluid"  world  threat.  In  that  past  decade, 
the  missile  it  was  initially  designed  to  destroy  has  become  faster,  smaller,  smarter 
and  more  maneuverable.  Furthermore,  the  U.  S.  fleet  has  also  made  a  transition 
from  operating  in  the  open,  "blue  water,"  sea  to  the  coastal,  "brown  water," 
environment.  With  this  change  in  operating  areas,  small  speed  boats,  which  were 
considered  an  insignificant  challenge  in  blue  water,  now  have  exposed  a  weakness 
in  surface  combatants.  Rather  than  building  a  new  weapon  system  to  deal  with 
this  threat,  the  U.  S.  Navy  has  decided  to  adapt  the  PHALANX  CIWS  into  an 
anti-patrol  boat  mode.  These  two  growing  threats,  with  their  decreasing  radar 
cross-section,  and  increasing  ability  to  pop-up  with  little  notice,  requires  the 
PHALANX  to  move  beyond  its  original  design  envelope.  The  primary  area  of 
concern  with  the  PHALANX  gun  is  to  tighten  its  dispersion  pattern  to  put  more 
rounds  on  target,  sooner  and  hence  destroy  the  target  at  a  greater  range  from  the 
ship.  This  is  the  focus  of  this  study. 

Before  discussing  the  nature  of  the  problem,  it  is  important  to  have  an  overall 
understanding  of  the  gun  system.  Figure  2-1  shows  a  component  breakdown  of  the 
PHALANX  weapon  system.  Starting  at  the  top,  the  dome  structure  contains  an  air 
search  radar  which  provides  the  initial  threat  detection.  It  has  a  maximum 
effective  range  for  detection  greater  than  4500  m  (5000  yds).  It  provides  360 
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degrees,  horizon  to  azimuth  search  coverage.  Once  a  target  has  been  identified  by 
the  search  radar  and  the  computer  determines  that  it  meets  the  criteria  of  a  threat 
to  the  ship,  the  gxm  moimt  slews  to  the  target  bearing  and  transfers  tracking  of  it 
to  the  fire  control  radar.  This  radar  is  mounted  below  the  search  radar,  on  the 
side  of  the  cylindrical  structure,  and  is  physically  boresighted  to  the  gun  boresight. 
It  is  a  high  resolution  radar  whose  purpose  is  to  track  the  target  as  well  as  the 
outgoing  rounds.  The  fire  control  computer  sends  the  necessary  servo  information 
to  keep  this  tracking  radar,  and  hence  the  gun,  locked  on  target  [Ref.  2]. 


Figure  2-1.  Subsystems  of  the  PHALANX  Close-in  Weapon  System  [Ref.  2]. 


Below  this  radar  is  the  heart  of  the  weapon  system.  It  is  the  six-barrel 
M61A1  gatling  machine  gun  as  shown  in  Figure  2-2.  It  fires  20  mm  (0.79  in) 


6 


tungsten  (formerly  depleted  uranium)  penetrators  at  a  velocity  of  1100  m/s  (3600 
fVsec).  An  example  of  a  projectile  is  shown  in  Figure  2-3.  The  penetrators  are 
designed  to  mostly  maintain  their  shapw  as  they  hit  and  pass  through  the  target. 
In  this  way,  the  penetrators  can  transfer  their  high  momentum  to  the  target  in  a 
concentrated  manner,  rather  than  fragment  and  perhaps  partially  reflect  off  the 
target.  In  experimental  studies  conducted  at  the  Naval  Surface  Warfare  Center  in 
Dahlgren,  Virginia,  it  takes  on  average  six  hits  to  destroy  an  incoming  missile. 
Tests  are  currently  being  conducted  against  patrol  boats,  but  no  data  is  available 
yet. 


Figure  2-2.  M61A1  PHALANX  gatling  machine  gun  [Ref  3.]. 
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Figure  2-3.  Mk  149  PHALANX  round  [Ref.  2]. 

The  gun  has  a  selectable  firing  rate  of  either  3000  shots  per  minute  or  4500 
shots  per  minute.  A  typical  missile  engagement  can  last  from  three  to  five  seconds 
so  that  the  higher  firing  rate  is  used  more  often.  This  means  that  somewhere 
between  225  and  375  bullets  are  fired  at  the  missile.  For  engagements  longer  than 
five  seconds,  the  barrels  begin  to  deteriorate  significantly  due  to  the  high  barrel 
temperature.  Therefore,  this  type  of  firing  burst  is  strongly  discouraged.  The 
ammunition  magazine  is  suspended  below  the  gun  and  has  been  upgraded  to  hold 
1500  rounds  allowing  for  roughly  three  engagements.  Only  one  barrel  fires  at  a 
time  near  the  top  of  the  barrel  group  while  the  spent  casing  extraction  and  breech 
loading  occurs  during  the  rest  of  the  rotation  cycle. 
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The  effective  range  of  the  gun  is  about  2000  m  (2187  yds).  However,  the  gun, 
after  calculating  the  target  speed,  begins  firing  when  the  target  is  about  5  seconds 
&om  the  firing  cutoff  point  at  100  m  (328  ft)  from  the  ship.  If  the  target  is  moving 
fast,  the  gun  will  fire  so  that  the  rounds  will  begin  intercepting  the  target  at  2000 
m.  If  the  target  is  moving  slow,  the  gim  will  wait  until  the  target  is  5  seconds  from 
the  ship  even  though  this  range  could  be  much  less  than  the  maximum  2000  m 
range.  The  firing  algorithm  will  be  different  for  the  anti-patrol  boat  surface  mode. 

Beneath  the  gun  is  a  cubical  structure  which  houses  the  computer,  radar 
electronics,  and  associated  cooling  systems.  This  is  the  brain  of  the  system.  It 
determines  the  threat,  time  to  impact,  where  to  point  the  gun,  when  to  start  firing, 
how  to  move  the  gun  based  on  observing  both  the  target  and  outgoing  rounds,  and 
when  to  cease  the  engagement.  It  is  a  closed-loop  spotting  system  in  that  the 
computer  receives  feedback  on  how  close  the  roiinds  are  coming  to  the  target  and 
calculates  the  necessary  gun  movements  to  ensure  intercept.  Figure  2-4  illustrates 
the  concept  of  close-loop  spotting.  The  fire  control  computer  calculates  aim  error  by 
measuring  the  projectile  miss  angle  at  the  target  and  subsequently  corrects  the 
gun’s  aimpoint.  When  the  fire  control  computer  determines  that  the  target  is 
destroyed,  that  is,  when  the  target’s  speed  or  radar  cross-section  drops  below  a 
certain  threshold,  then  the  firing  stops  and  the  gun  slews  back  to  its  horizontal 
position  and  resumes  the  search  mode.  For  ships  with  two  or  more  mounts,  they 
are  able  to  talk  to  each  other  in  order  to  prevent  engaging  the  same  target.  [Ref. 
2] 

PHALANX  defines  an  air  target  as  something  heading  towards  the  ship,  or 
something  that  could  be  heading  towards  the  ship  after  a  9  G  turn.  This  precludes 
the  possibility  that  the  gun  will  shoot  down  the  ship’s  own  outbound  missiles.  The 
speed  of  the  target  has  to  be  above  a  certain  limit,  so  that  friendly  helicopters,  and 
such,  are  not  accidentally  engaged.  The  PHALANX  gun  is  designed  to  protect  the 
ship  on  which  it  is  mounted,  so  it  will  not  attempt  to  defend  other  ships. 
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Figure  2*4.  Geometry  of  PHALANX  missile  engagement  [Ref.  4]. 

B.  PHALANX  ROUND  DISPERSION 

Based  on  dispersion  data  provided  during  a  data  gathering  trip  to  Hughes 
Missile  Systems  Division,  Pomona,  California,  it  was  learned  that  the  PHALANX 
weapon  system  knows  where  the  target  is,  but  has  had  a  hard  time  hitting  it 
effectively  as  evidenced  by  the  dispersion  pattern  shown  in  Figure  2-5.  Here,  one 
can  see  the  six  individual  barrel  dispersion  patterns,  which  are  generally  good 
individually.  However,  the  target  is  the  object  at  the  center  of  this  particular 
dispersion  pattern  where  virtually  no  rounds  have  fallen.  The  fire  control  system 
is  doing  its  job  by  keeping  the  target  at  the  center  of  the  gun’s  dispersion  pattern. 
It  appears  that  much  of  the  dispers'‘on  can  be  associated  with  the  six  barrel-gun 
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system.  To  elaborate  on  dispersion.  Figure  2-6  shows  the  geometry  used  to  define 
dispersion.  The  probability  of  a  penetrator  passing  through  a  circle  of  radius  r  and 
standard  deviation  o  is  given  by  the  distribution, 


p(r)  = 


-  2o* 
e 

2na^ 


(2.1) 


It  is  assumed  that  all  the  penetrators  pass  through  a  circle  of  infinite  radius,  so 
that, 

2ii 

rp(r)dnci8  =  1  (2.2) 

The  probability  that  a  penetrator  will  pass  within  a  circle  of  radius  one  o,  P^g,  can 
be  determined  by  integrating  equation  (2.1)  in  polar  coordinates  from  r  =0  to  r  =  o. 


P<a  =  l^pirydrdB 


^0  2ro* 


rdr2ji 


(2.3) 


Simplifying  and  solving  the  equation  yields, 


=  +  1  =  0.39  . 


Therefore,  39  percent  of  the  penetrators  fall  withing  one  a  while  61  percent  fall 
outside.  In  the  dispersion  example  of  Figure  2-6,  39  percent  of  the  penetrators  of  a 
PHALANX  gun  will  pass  within  a  circle  of  one  o  =  1.2  milliradians.  Part  of  the 
dispersion  ('1  milliradian)  can  be  attributed  to  spin  induced  yaw  after  the  biillet 
leaves  the  barrel.  But,  a  significant  part  of  the  dispersion  is  directly  attributed  to 
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barrel  tip  motion.  The  observed  2.4  milliradians  of  dispersion  corresponds  to  3.6 
mm  of  barrel  tip  displacement  of  a  1.52  m  (59.8  in)  long  barrel. 


.TOOPBeSENTEDAREAAT  AZIMUTH  (MRAO> 

(000  R  WITH  200  PT  OFF  RANGE 


Figure  2*5.  Experimental  dispersion  pattern  of  the  PHALANX  gun.  Note 

that  the  target  drone  is  at  the  center  while  the  individual  barrel 
dispersion  patterns  fall  all  around  it  [Ref  3]. 

Other  dispersion  patterns  shown  in  Figure  2*7  show  a  less  pronounced 
doughnut-shaped  pattern  than  in  Figure  2-5,  but  it  is  still  generally  visible,  with 
the  target  remaining  unscathed  in  the  center.  The  dispersion  patterns  were 
measured  using  an  Oehler  and  battenboard  system  at  a  test  facility  in  China  Lake 
as  shown  in  Figure  2-8.  An  Oehler  system  consists  of  a  set  of  four  microphones 
which  will  listen  for  the  acoustic  shock  wave  of  the  parsing  supersonic  projectile. 
The  arrival  time  of  each  microphone  is  compared  with  the  other  microphones  and  a 
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Figure  2-6.  Illustration  of  radial  dispersion  as  defined  for  the  PHALA>TX  gun 
[Ref.  3]. 

bullet  position  is  calculated.  The  battenboard  follows  the  Oehler  system  and  is 
used  to  record  the  roimd  pattern  and  to  distinguish  between  the  pusher  plugs, 
which  drive  the  ro\ind  through  the  barrel  (see  Figure  2-3),  and  the  rounds 
themselves.  This  observation  set-up  was  designed  to  test  the  gun  dispersion  and 
not  that  of  the  ammunition. 

The  PHALANX  dispersion  pattern  measurements  also  show  that  the  six- 
barrel  system  gives  a  dispersion  pattern  greater  than  that  of  a  single  barrel  fired 
separately.  Thus,  the  structural  dynamics  of  the  six-barrel  system  appears  to 
warrant  an  in-depth  analysis.  According  to  the  Hu|^es  data,  the  PHALANX  gun 
system  has  a  typical  dispersion  range  of  1.2  to  1.5  milliradians  and  a  maximum  of 
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Figure  2-7.  A  less  pronounced  circular  PHALANX  dispersion  pattern  [Ref  3]. 

1.6  milliradians.  The  U.  S.  Navy  would  like  to  achieve  an  average  dispersion  of 
less  than  1.2  milliradians.  In  order  to  achieve  this  goal,  the  dynamics  of  the  gun 
system  will  have  to  be  better  understood  before  any  type  of  fix  can  be  fully 
implemented.  This  is  the  premise  behind  this  study.  What  effect  the  observed 
dispersion  has  on  target  destruction  is  the  subject  of  the  next  section  and  was  the 
motivation  for  imdertaking  this  project. 
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Figure  2-8.  Experimental  set-up  for  observing  the  PHALANX  dispersion 
patterns  [Ref.  3]. 

C.  THE  BALLISTIC  DEBRIS  PROBLEM 

Although  the  maximum  effective  range  of  the  PHALANX  g\m  is  2000  m,  field 
and  computer  simulations  have  shown  that  the  destruction  of  an  incoming  missile 
does  not  occur,  if  at  all,  imtil  after  the  missile  is  within  1(X)0  m  (3281  ft)  of  the 
ship.  On  the  surface  this  does  not  appear  to  be  a  problem.  As  part  of  a  class 
computer  simulation  project,  the  resulting  debris  firom  the  destroyed  missile 
continues  on  a  ballistic  tr^ectory  toward  the  ship  at  roughly  the  speed  of  the 
missile.  The  firagment  sizes,  which  can  be  on  the  order  of  a  few  to  a  hundred 
pounds,  coupled  with  their  high  velocity,  which  is  typically  from  1000  m/s  to  1500 
m/s,  would  cause  a  great  deal  of  damage  to  an  unarmored  superstructure  which 


15 


usually  houses  sensitive  electronics,  antennas  and  not  to  mention  sailors.  Metallic 
particles  at  1000  m/s  would  easily  punch  holes  through  several  bulkheads  before 
being  stopped.  Although  the  PHALANX  gun  "worked"  in  destrojdng  the  missile, 
the  ship  would  still  suffer  tremendous  damage.  A  solution  to  this  problem,  would 
be  to  obviously  destroy  the  missile  at  a  greater  distance  from  the  ship.  In  order  to 
achieve  this  goal  the  dispersion  pattern  would  have  to  be  tightened  so  more  hits 
occurred  at  near  the  beginning  of  the  PHALANX  engagement  rather  than  at  the 
end. 

To  clarify  these  points  further,  a  computer  simulation  is  used  to  describe 
fragment  trajectories.  In  this  model,  a  missile  is  assumed  to  be  destroyed  after  one 
or  several  of  the  PHALANX  rounds  detonates  the  warhead  of  the  missile.  The 
range  at  which  this  detonation  occurs  is  varied  so  that  the  number  of  fragments 
that  hit  the  ship  can  be  determined.  The  fragment  size,  shape  and  mass  is 
approximated  by  distributing  the  fragment  terminal  velocities  by  the  exponential 
distribution. 


where  is  a  fragment  terminal  velocity,  and  is  the  characteristic  terminal 
velocity  which  is  derived  frrom  typical  fragment  attributes.  The  terminal  velocity  is 
calculated  by  the  expression 


(2.5) 


where  m  is  the  mass  of  the  fragment,  g  is  the  acceleration  of  gravity,  C  is  the 
coefEcient  of  air  drag,  p  is  the  density  of  air  and  A  is  the  debris  fragment  surface 
area.  Additional  velocity  components  caused  by  the  force  of  the  explosion  are 
assigned  to  the  fragments  along  each  of  the  three  coordinate  axes.  These  velocities 
are  determined  from  a  Gaussian  distribution  with  a  standard  deviation  of  0^,  in 
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each  direction.  The  value  of  =  10  m/s  was  provided  by  the  Naval  Surface 
Warfare  Center,  Dahlgren.  When  the  values  of  the  variables  described  above  are 
entered  into  the  simulation  program,  the  results  appear  like  that  of  Figure  2-9 
This  figure  shows  a  typical  best  case  scenario  where  the  missile  blew  up  at  the 
maximum  PHALANX  range  of  2000  m.  It  is  divided  up  into  a  side  view,  top  view, 
and  an  end  view  showing  how  the  missile  disintegrates  as  viewed  from  the  ship. 
The  missile  is  destroyed  at  the  left  most  edge  in  the  top  and  side  views.  The 
"ship"  is  represented  as  the  right  hand  edge  of  these  views.  The  black  dots  are  the 
missile  debris  fragments  and  the  black  lines  represent  their  tr^ectories.  The 
"ocean"  is  at  the  bottom  edge  of  each  view.  It  is  easy  to  record  ship  hits  as  those 
that  reach  the  right  hand  edge.  In  this  particular  case,  the  missile  was  traveling 
at  a  fast  speed  of  =  1000  m/s  (3281  fl/s)  and  a  height  of  10  m  above  the  water 
when  it  exploded.  NF  stands  for  the  number  of  fragments.  Here  6  missile 
firagments  hit  the  ship  with  a  velocity  of  around  100  m/s  (328  fVs). 

In  Figxire  2-10,  the  missile  has  the  same  characteristics,  but  now  does  not 
explode  until  it  is  intercepted  by  the  last  PHALANX  bullet  at  a  range  R  =  500  m 
(1641  ft)  from  the  ship.  The  resulting  debris  hits  the  ship  in  a  "shotgun  blast" 
pattern  where  18  of  the  20  fragments  reach  the  ship.  This  latter  case  is  the  more 
likely  of  the  two  scenarios  for  a  nonmaneuvering,  medium  speed  missile  [Ref  5]. 

If  the  missile  was  more  sophisticated  in  that  it  "jinked,"  or  performed  random 
maneuvers,  or  if  it  had  a  higher  speed,  then  the  intercept  range  is  even  more 
reduced.  Therefore,  it  is  readily  apparent  that  significant  debris  from  an  effectively 
engaged  missile  is  going  to  hit  the  defending  ship  in  virtually  every  case.  Such 
debris  would  certainly  not  destroy  the  ship,  but  might  seriously  impair  her  ability 
to  fight. 

In  other  simulations,  performed  by  Professor  W.  B.  Colson,  at  the  Naval 
Postgraduate  School,  the  effectiveness  of  the  gun  in  destroying  an  incoming  missile 
is  tested.  In  this  simiilation,  the  round  dispersion  is  assumed  to  have  a  Gaussian 
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Figure  2-9.  TMgectories  of  fragments  of  a  missile  destroyed  by  a  PHALANX 
gun  at  a  range  R  =  2000  m.  Notice  6  of  the  20  missile  fragments 
hit  the  ship  which  is  at  the  right  hand  edge  of  each  view. 

distribution  with  a  standard  deviation  using  a  typical  PHALANX  dispersion  value 
of  50  =  2  milliradians.  Figure  2-11  shows  the  results  of  a  PHALANX  engagement 
with  an  incoming  missile  which  has  a  speed  of  =  600  m/s  (1969  fVs)  and  a 
physical  cross-section  of  0.20  m  (7.9  in).  In  thiR  scenario,  the  missile  speed  is  not 
unusual,  but  it  does  begin  a  terminal  5  G  maneuver  at  a  range  of  roughly  2000  m 
(6562  ft).  The  PHALANX  gun  destroys  the  missile  at  a  range  of  roughly  200  m 
(656  ft)  after  it  was  only  hit  three  times.  The  number  of  hits  required  to  destroy 
the  missile  is  chosen  randomly  and  is  better  than  the  usual  six.  The  effect  of  the 
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Figure  2-10.  An  incoming  missile  is  destroyed  by  the  PHALANX  gun  at  a 
range  of  500m.  Most  of  its  debris  hits  the  ship. 

maneuver  was  to  allow  the  missile  to  penetrate  to  within  200  m  of  the  ship  when  it 
exploded  and  sends  almost  all  of  its  debris  at  the  ship. 

Another  engagement  is  shown  in  Figure  2*12  where  the  missile  does  not  have 
a  terminal  maneuver,  but  has  a  higher  speed  of  800  m/s  (2625  ft/s).  In  this  case, 
the  missile  is  also  destroyed,  but  like  before,  it  is  not  destroyed  far  enough  away 
from  the  ship  so  all  of  its  debris  strikes  the  ship.  [Ref  5] 
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Figure  2-11.  A  terminal  missile  maneuver  defeats  the  PHALANX  gun  [Ref.  5], 

These  studies  provided  the  motivation  for  trying  to  improve  the  dispersion  of 
the  PHALANX  gun.  In  order  to  make  such  improvements,  the  gun  dynamics  have 
to  be  first  understood.  This  study  attempts  to  provide  a  solid  foundation  upon 
which  to  study  the  gun.  In  an  era  of  reduced  defense  spending,  but  not  necessarily 
one  of  reduced  potential  for  hostilities,  the  PHALANX  gun  is  going  to  remain  a 
ship’s  last  layer  of  defense  into  the  next  century. 
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Figure  2-12.  A  high  speed,  nonmaneuvering  missile  defeats  the  PHALANX 
gun  [Ref.  5]. 
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III.  FINITE  ELEMENT  THEORY 


A.  THE  DECISION  TO  USE  FINITE  ELEMENT  MODELING 

In  order  to  investigate  the  mechanical  dynamics  behind  the  dispersion  of  the 
PHALANX  six-barrel  system  and  to  study  the  various  design  factors  which 
influence  projectile  trajectory  and  accuracy,  it  was  decided  that  a  finite  element 
computer  model  would  be  the  best  method  for  gaining  insi^^t  into  this  complex 
system.  Many  aspects  of  the  gun  dynamics  would  have  been  difficult  if  not 
impossible  to  understand  by  studying  the  gun  itself  during  its  operation  due  to  its 
hiyh  rate  of  fire  and  extremely  high  barrel  temperatures  as  well  as  with  problems 
of  on-sight  availability.  Even  high  speed  film  does  not  provide  as  much  benefit  or 
nearly  the  flexibility  as  a  computer  model.  A  finite  element  analysis  can  provide  a 
thorough  model  of  the  gun  where  various  driving  forces  can  be  applied  and  the 
gun’s  response,  such  as  barrel  tip  motion,  calculated.  This  also  allows  individual 
components  to  be  studied  in  detail  to  understand  their  affect  on  the  overall 
dispersion  pattern.  Furthermore,  this  finite  element  computer  model  is  adaptable 
enough  to  support  redesigning  the  PHALANX  barrels  and  perhaps  the  thrust 
bearing,  as  well  as  further  modifications,  such  as  a  change  in  projectile  size  or  the 
addition  of  a  gun  barrel  restraint.  Finally,  a  computer  code  provides  an  economical 
solution  by  leapfi:x)gging  past  design  changes  made  possible  by  computer  simulation 
instead  of  expensive  hardware  tests. 

The  finite  element  computer  code  called  I-DEAS  (Integrated  Design 
Engineering  Analysis)  developed  by  SDRC  (Structural  D}mamics  Research 
Corporation)  was  chosen  to  model  the  PHALANX  gun  due  to  its  existing  on-site 
license  with  the  Naval  Postgraduate  School  and  its  "completeness."  That  is,  it 
precluded  the  need  for  other  software  packages  allowing  for  a  much  smoother 
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transition  between  model  construction  and  testing.  I-DEAS  is  made  up  of  the 
following  set  of  modules:  Solid  Modeling,  Drafting.  Finite  Element  Modeling  and 
Analysis,  System  Dynamics,  Test  Data  Analysis,  and  Manufactiuing.  This  study 
uses  a  two-step  approach.  The  first  step  uses  the  Finite  Element  Modeling  and 
Analysis  Module  to  develop  a  finite  element  model  of  the  gun.  This  module  allows 
normal  modes,  stresses,  and  static  defiections  to  be  calculated  after  the  structure’s 
material  properties  and  loads  are  defined.  In  the  second  step,  the  completed  finite 
element  model  is  transferred  to  the  System  Dynamics  Analysis  module  where  an 
in-depth  forced  response  analysis  is  then  performed. 

The  primary  purpose  behind  the  finite  element  modeling  of  the  PHALANX 
gun  is  to  better  understand  some  of  the  causes  behind  the  apparent  dispersion 
pattern.  Specifically,  barrel  tip  displacement  pla3rs  a  critical  role  in  the  dispersion. 
One  of  the  best  ways  to  predict  the  motion  of  a  system  while  being  excited  is  to 
first  solve  for  its  normal  modes.  Normal  modes  of  vibration  have  the  iiseful 
property  of  being  able  to  describe  any  complex  motion  (due  to  arbitrary  initial 
conditions)  as  a  superposition  of  normal  mode  amplitudes.  So  what  are  normal 
modes  of  vibration?  Briefly,  every  structure,  with  its  inherent  mass  and  stiflhess, 
has  its  own  distinctive  set  of  normal  modes.  A  structure  vibrating  at  one  of  its 
natural  fii^quendes  does  so  with  a  characteristic  pattern  of  amplitude  distribution 
called  a  normal  mode  of  vibration.  A  more  mathematical  explanation  will  be 
provided  as  the  theory  behind  finite  elements  unfolds  in  this  chapter.  After 
modeling  a  structure  as  a  series  of  finite  elements  and  connecting  points,  or  nodes, 
the  computer  code  solves  for  these  normal  modes  of  vibration.  Then,  this  model 
with  its  normal  modes  is  irm^iiferred  to  the  Systems  Dynamics  Analysis  module 
where  excitation  forces  and  damping  can  be  applied.  The  response,  such  as  barrel 
tip  displacement,  can  then  be  calculated.  A  fourier  transform  of  the  time  response, 
which  performs  a  transient  integration,  estimates  which  normal  modes  are  excited 
and  their  relative  degrees  of  excitation.  This  same  information  can  be  revealed 
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using  a  frequency  response  function  which  is  the  method  presented  in  Chapter  VI. 
This  latter  method  has  the  added  benefit  of  not  requiring  a  fourier  transform.  The 
resulting  dispersion  could  then  be  modified,  for  example,  by  tuning  these  normal 
modes  to  different  frequencies. 

The  finite  element  method  provides  a  discrete  representation  of  a  dynamical 
system  [Ref.  6].  It  involves  the  systematic  solution  of  partial  differential  equations 
allowing  the  structural  analysis  of  the  complex  systems  to  be  easily  performed  by 
digital  computers.  Developed  in  the  late  1950’s,  it  has  been  made  more  accessible 
with  the  advent  of  powerful  desktop  computers.  It  is  especially  iiseful  for 
structures  that  contain  a  large  number  of  degrees  of  fii^dom,  where  a  variety  of 
materials  are  used  or  where  discontinuities  such  as  thickness  variations,  holes,  or 
curved  beams  are  present  [Ref  7].  These  types  of  structures  make  exact  solutions 
virtually  impossible. 

The  finite  element  method  breaks  down  a  structure  into  simple,  discrete 
pieces  such  as  beams,  plates,  or  bricks,  each  of  which  is  responsible  for  carrying 
some  portion  of  the  load  applied  to  the  structure.  The  equations  describing  the 
motion  or  elasticity  of  these  discrete  elements  can  be  solved  individually  with 
boundary  conditions  matched  at  their  nodal  points  to  get  the  overall  dynamic 
response  of  the  structure.  The  calculated  solution  is  not  exact,  but,  with  good 
modeling  practice  the  error  can  be  reduced  to  tolerable  levels.  Moreover,  there  is 
no  other  way  to  solve  the  problem.  For  displacement  based  finite  element 
modeling,  the  error  decreases  monotonically.  The  finite  elements  and  nodes  are 
assembled  to  correspond  as  close  as  possible  to  the  geometry  of  the  actual 
structure  being  modeled.  More  elements  are  assembled  in  areas  where  the 
greatest  stress  gradient  or  curvature  is  expected  or  where  the  greatest  accuracy  of 
solution  is  needed.  This  often  leads  to  a  trade-off  between  computational  speed 
and  computer  storage  on  the  one  hand  and  the  necessary  accuracy  on  the  other. 
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B.  THE  DISCRETIZATION  OF  A  SIMPLE  ROD 


As  an  illustration  of  the  finite  element  method,  consider  a  single  rod  element 
of  length  L  as  shown  in  Figure  3.1  below.  It  has  a  node  at  each  end  point  whose 
coordinates  in  the  element’s  local  coordinate  system  (x  ^ )  are  x  =  0  for  node  1  and 
X  sL  for  node  2.  This  element  is  situated  in  a  global  coordinate  system  (x,z) 
which  would  be  used  to  describe  a  structure  of  elements. 


I  M  I 


Figiire  3-1.  A  single  rod  element. 

To  structurally  represent  the  rod’s  continuous  distribution  of  mass,  it  will  be 
broken  down  into  an  infinite  number  of  discretized  mass  units  dm .  These  mass 
units  are  then  idealized  by  tracking  their  positions  to  the  coordinates  at  the  end  of 
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the  rod,  specifically  its  nodes.  Finally,  an  axial  force  is  applied  to  the  element 
along  the  rod’s  local  x  •  axis  which  will  cause  an  axial  deformation. 

In  order  to  investigate  the  forces  and  displacements  of  this  rod  element  it  will 
be  beneficial  to  use  a  displacement  field  variable  (  u  =u{.x^))  to  describe  the 
positions  of  the  nodes  at  any  instant  in  time.  This  is  important  since  it  is  at  these 
node  points  that  the  rod  element  is  coupled  to  the  rest  of  the  structure.  Any  load 
information  transmitted  firom  the  rest  of  the  structure  is  applied  to  the  rod  element 
at  its  nodes.  Therefore,  the  response  anywhere  along  the  rod  can  be  described  in 
terms  of  the  response  at  the  element’s  endpoints,  that  is,  its  nodal  degrees  of 
fi’eedom,  which  are  finite.  This  is  a  finite  discretization  of  an  infinite  element 
system.  An  exact  solution  of  the  displacement  field  is  difficult  considering  that  the 
force  and  moment  distribution  as  well  aa  the  botindary  conditions  at  the  ends  of 
the  rod  are  not  known.  However,  an  approximation  can  be  made  by  considering  a 
finite  number  of  degrees  of  fi:eedom  for  the  displacement  fields  of  each  element 
[Ref.  7].  The  coordinates  for  nodes  1  and  2  at  any  instant  in  time  are  respectively, 

ui^um)  (3.1) 

and 

U2-u{Ljt)  .  (3.2) 

Since  the  load  is  an  axial  force,  only  displacements  along  the  loc*!  x  •  axis  are 
considered,  so  the  rod  deformation  can  be  described  by  these  two  displacement 
degrees  of  fireedom,  u  x  and  u  2- 

The  displacement  field  is  a  linear  motion  in  x  since  there  are  two  axial 
degrees  of  fii^edom  [Ref.  7].  It  is  written  as 

u(x,t )  =  Cl  +  C2X  ,  (3.3) 

where  Ci  and  C2  are  vinknown  constants.  Substituting  equation  (3.3)  into  equations 
(3.1)  and  (3.2)  yields 
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«i  =  Ci 


(3.4) 


U2  =  Ci+C2^  ,  (3.5) 

which  forms  a  system  of  simultaneous  equations 


1  0 

Cl 

1  L 

.^2. 

.“2 

(3.6) 


where 


1  0 

-1 

1 

L  0 

1  L 

”  L 

-1  1 

Solving  for  the  Cj  coefficients  gives 

Ci  =  Ui 


and 


(3.7) 


(3.8) 


Cj  = 


(— tt  1  +  ^2) 

L 


Substituting  equations  (3.8)  and  (3.9)  into  equation  (3.3) 


u(x^)  =  ui 


L  L 


(3.9) 


,  *1 

X 

‘-T 

Ui  + 

L 

»  4 

=  Oi(x  )u  i« )  +  a2(x  )u  2(t )  .  (3.10) 


where 


and 


(3.11) 


)=j-  , 


(3.12) 


or  in  matrix  notation 
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where 


u  (x 


.1: 


“1 

“2 


(3.13) 


a 


ajCx)!' 


(3.14) 


The  response  anywhere  along  the  rod  can  now  he  described  in  terms  of  the  nodal 
coordinates  and  the  interpolation  functions  or  shape  functions  ot^ix).  These 
functions  interpolate  the  element  displacement  between  the  nodes  according  to  the 
boundary  conditions  and  movement  of  the  nodes.  They  also  form  a  basis  set  where 
any  motion  of  the  rod  can  be  described  as  a  linear  combination  of  these  shape 
functions. 

The  boimdary  conditions  defined  earlier  (see  Figure  3-1)  for  the  element 
impose  the  following  requirements  on  the  shape  functions,  Oj : 

ai(0)  =  1  aa(0)  =  0  (3.15) 

ai(L)*0  aa(L)=l  ,  (3.16) 

where  each  shape  function  represents  the  element  shape  when  a  unit  displacement 
is  applied.  Looking  at  Figure  3*2,  it  is  apparent  that  the  shape  functions  can  now 
be  written  explicitly  as 

ai(x )  =  1-L 


a^x )  *  ^ 

The  shape  fimctions  do  not  give  an  exact  solution  of  the  displacement  field 
between  the  nodes,  but  as  the  niunber  of  elements  are  increased  along  the  length 
of  the  rod  a  better  solution  is  obtained.  Referring  to  Figure  3-3,  an  exact  solution 
of  the  displacement  field  can  be  written  as 
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a^L )  =  1 


ai(x )  =  1  -  L 

ai(0)  =  1 


0  X 


Figure  3-2.  Shape  functions  . 


(3.17) 


where  x’  and  y’  are  the  coordinates  of  a  point  P  along  the  length  of  the  rod  and  the 


quantities 


and 


are  the  basis  vectors.  Therefore,  any  point  along  the  length 


of  the  continuous  rod  can  be  described  in  terms  of  the  basis  vectors.  However,  it  is 
the  quantities  x’  and  y’  which  are  unknown  and  must  be  determined  according  to 
the  motion  of  its  endpoints.  This  is  a  linear  assumption  that  is  only  valid  if  the 
transverse  displacement  is  small,  that  is,  x/Lcl. 
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Figure  3-3.  Diagram  of  a  point  P  along  a  deformed  rod. 

C.  MATRIX  FORM  OF  EQUILIBRIUM  EQUATIONS 

Since  each  normal  mode  of  vibration  is  made  up  of  the  periodic  exchange  of 
kinetic  and  potential  energies  within  the  structure,  it  is  necessary  to  determine  the 
kinetic  and  potential  energies  for  each  element  [Ref.  8].  Specifically,  the  mass  and 
stiffiiess  matrices  of  the  element  can  be  deduced  after  applying  Lagrange’s 
equation  to  the  internal  kinetic  and  potential  energies.  In  the  case  of  the  rod,  its 
kinetic  energy  T  and  strain  energy  U  are  respectively 

T  =  ^j^[u(x^)]^dx  (3.18) 

TT  _  EA.  f  ^  du{xjt) 

"  2  dx 


(3.19) 


where  y  is  the  density  per  unit  length  and  E  is  Young’s  moduliis  of  elasticity  and  it 
is  assumed  that  the  rod  is  uniform  with  a  constant  cross-sectional  area  A .  The 
strain  energy  occurs  due  to  any  axial  elongation  of  the  rod  where  du(x  ,t)  i  dx  is  its 
corresponding  strain  rate  of  change. 

Differentiating  equation  (3.10)  with  respect  to  time  yields 

u  =  ai(x)ui(/)  +  a2(*)ji2«)  (3.20) 

and  with  respect  to  z  gives 

4^  =  a'i(x)ui(0  +  a'2(i)tt2(«)  •  (3.21) 

dx 

Now  substituting  these  expressions  back  into  equations  (3.18)  and  (3.19)  defines 
the  rod’s  internal  kinetic  and  potential  energies  respectively  as 

V 

T  =  (oi(x  )u  1  +  02(x  liiate  (3.22) 

t;  =  M  J^V'i(x  )u  X  +  aVx  )u  ifdx  .  (3.23) 

The  equations  of  motion  can  be  derived  firom  T  and  U  by  applying  Lagrange’s 
equation.  It  is  defined  as 

=  0  ,  (i  =  1,2,3, ...n)  .  (3.24) 

dt  dq.  9(7, 

where  L  is  the  Lagrangian  of  the  system.  For  conservative  systems,  where  the 
sum  of  potential  and  kinetic  energies  remains  constant,  and  mechanisms  of  energy 
loss  such  as  internal  firiction  are  ignored,  the  Lagrangian  is  the  difference  between 
kinetic  and  potential  energies: 

L  =  T  -U  .  (3.25) 

It  is  important  to  note  that  damping  as  a  form  of  energy  loss  can  be  added  into  the 
equations  of  motion  later  (see  Chapter  VI).  The  generalized  coordinates  Qi 
(i=l,2,3,...n )  describe  the  configuration  of  the  system,  with  n  degrees  of  freedom 
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and  hence  n  equations  [Kef.  9].  In  the  case  of  the  rod  element,  the  generalized 
coordinates  are 


=  Ui  Q2  -  ^2  > 


so  that 


(3.26) 


Now  applying  Lagrange’s  equation 

[ai(i  )u  1  +  02(1  )ii  2]oii(i  )dx  (3.27) 

=  EAj  [a'i(x  )u  1  +  aVx  )u  2)0' i(z  )dx  (3.28) 

au  \ 

Therefore,  the  element  mass  matrix  coefficients  can  be  written  as 

niij  =  Oi  (x  )ajix  )dx  ,  (3.29) 

and  the  stiffiaess  coefficients  matrix  as 

kij  =EAj^a'i(x)a'j(x)dx  .  (3.30) 

So  the  equations  of  motion  can  be  written  in  the  familiar  Newtonian  mechanical 
forms 


>l 

*11 

*12 

Ml 

F2 

• 

*21 

*22 

.“2 

•  ■ 

^1 

.-^2 

= 

mil  mj2 

m  21  w  22 

Ml 

“2 

(3.31) 


(3.32) 


It  is  important  to  note  that  equations  (3.29)  and  (3.30)  have  the  general 
quadratic  form  with  n  variables  qi  of  the  type 


(3.33) 


where  the  aij  (i.e.  y  and  EA)  are  coeffidenta  independent  of  [Ref.  9].  This  form 
is  diagonal  if  it  contains  terms  with  squares  of  Qi  namely  anq;,  022^2  >  > 

no  cross  terms  like  a  i2<7 19  2  9]- 

However,  if  there  are  cross  terms  in  either  the  mass  or  stiilhess  matrices, 
then  there  are  a  set  of  n  coupled  differential  equations  which  describe  the  motion 
of  the  system.  In  the  physical  sense,  this  means  that  the  motion  of  one  node 
influences  the  motion  of  the  other  through  these  coupling  terms.  They  are  linearly 
independent.  So  the  problem  now  is  to  simultaneously  diagonalize  the  mass  and 
stifBaess  matrices  if  they  are  not  already  diagonalized.  It  is  possible  to  exdte  a 
normal  mode  of  one  node  withou  ,  exdting  any  other  modes.  The  set  of  vectors  that 
diagonalize  the  modal  matrix  are  called  "normal  modes."  This  diagonalization  is 
achieved  by  performing  a  linear  transformation  involving  the  orthogonal  modal 
matrix  0,  whose  columns  are  the  eigenvectors  of  the  matrix  K  with  respect  to  M  of 
the  eigenproblem 

iK  -  ©2^)0  =  0  .  (3.34) 

Once  the  matrices  are  diagonalized,  and  hence  uncoupled,  then  the  n  equations  of 
motion  are  easily  solved  for  less  than  n,  rather  than  solving  for  n  coupled 
(simultaneous)  differential  equations,  and  more  importantly,  each  uncoupled 
differential  equation  is  one  of  the  normal  modes.  This  process  is  called  a  transition 
to  normal  coordinates.  It  is  important  to  note  that  the  physical  nature  of  the 
problem  has  not  been  changed,  rather,  the  normal  coordinates  merely  capture  the 
physical  configuration  of  the  system  so  as  to  greatly  simplify  the  mathematical 
description  of  the  system. 
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D.  FREE  VIBRATION  IN  A  STRUCTURE  OF  ELEMENTS 


Now  that  the  equations  of  motion  of  a  single  rod  element  are  understood,  a 
whole  structure  can  now  be  composed  of  elements  like  these.  The  equation  of 
motion  for  an  assembly  of  element  matrices  is 

[Mli+[iini=0  .  (3.35) 


where  [M]  and  [K]  are  square  matrices  and  the  quantities  with  a  single  underline 
imply  a  column  matrix.  As  an  example,  if  the  original  continuous  rod  is 
represented  by  eight  connected  elements  instead  of  one,  as  shown  in  Figure  3-4, 
each  one  with  its  own  2x2  stiffness  matrix,  Ki ,  then  the  eqxiation  of  motion  would 
appear  as 


\k,]  .  .  .  ■ 

“1 

F 

.  [^23  .  . 

“2 

L  J 

•  *  •  • 

• 

.  .  [Xg] 

I 

(3.36) 


where  each  element  stif&ess  matrix,  X,-,  has  an  overlapping  stifiness  with  adjacent 
elements  at  their  common  degrees  of  freedom. 

The  general  method  for  solving  the  free  vibration  problem,  that  is,  finding  the 
eigensolution  of  (3.35)  is  to  use  a  solution  of  the  form 


x«)  =  Cie*“' 


where  x(t)  is  the  column  matrix 


xit)  = 


Xiit) 

x^t) 


(3.37) 


(3.38) 


describing  the  position  of  the  element  end  points  or  nodes.  The  solution  consists  of 
a  vector  spatial  part. 


$  = 


♦i 

♦2 


(3.39) 
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Figure  3*4.  Continuous  rod  divided  into  8  elements. 


scalar  temporal  part,  e‘“*  and  a  complex  constant  C.  This  form  of  solution  was 
used  since  the  mass  and  stiffness  associated  with  each  element  is  similar  to  a 
mass/spring  system  which  executes  simple  harmonic  motion. 

Taking  the  necessary  derivatives  of  the  solution  and  substituting  back  into 
the  equation  of  motion  yields  the  homogeneous  equation 

(-a)2[Af  1$  +  =  0  .  (3.40) 

If  C  =  0  then  there  would  be  no  motion,  so  that  possible  solutions  to  this  equation 
have  to  come  firom 


([Jf]-a)2[Af])$  =  0  , 


(3.41) 


which  is  an  eigenvalue  problem.  Part  of  the  process  in  finding  solutions  to 
eigenvalue  problems  is  to  take  its  determinant 

det([K]-<aHM])^0  ,  (3.42) 

which  is  called  the  "characteristic  equation."  This  is  the  necessary  and  sufficient 
condition  for  these  n  homogeneous  equations  in  n  unknowns  to  possess  a  non¬ 
trivial  solution.  The  values  of  which  render  ([/f]  -  (o*[Af ])"^  singular  are  called 
eigenvalues.  For  two  degrees  of  freedom,  there  will  be  two  roots  of  this  equation, 
coi  and  (Os.  The  nodes  can  vibrate  in  simple  harmonic  motion  with  only  these  two 
frequencies.  To  determine  the  actual  motion  of  the  system,  Oi  and  (O2  are 
substituted  back  into  equation  (3.37)  to  find  Xi(t)  and  X2(t)  for  each  frequency. 
Since  the  equation  of  motion  is  homogeneous,  the  coefficients  Ci  and  C2  are 
undetermined  which  means  the  system  «in  vibrate  with  an  arbitrary  amplitude  for 
either  characteristic  fi?equency.  The  relation  between  Ci  and  C2  is  essentially  a 
phase  relation.  [Ref.  9] 

The  next  step  is  to  apply  this  method  to  the  whole  system  of  elements  with  n 
degrees  of  fireedom.  Starting  with  the  eigenproblem  defined  earlier 

=  0  ,  (3.41) 

except  now  [K]  and  [M]  are  made  up  of  a  whole  series  of  two  by  two  matrices,  one 
for  each  element.  For  systems  with  many  degrees  of  fireedom  and  especially  with 
cross  terms  appearing  in  the  mass  and  stifEhess  matrices,  the  n  coupled 
differential  equations  describing  the  motion  now  become  difficult  to  solve. 
However,  a  relatively  easy  method  to  solve  such  a  complex  system  is  to 
simultaneously  diagonalize  the  matrices  [M]  and  [ff].  Then,  it  is  trivial  to  find  the 
eigenvalues  and  eigenvectors.  In  this  case,  the  diagonal  elements  are  now  the 
eigenvalues  and  the  eigenvectors  are  then  the  columns  of  the  matrix  used  to 
diagonalize  [M]  and  [iT]. 
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The  method  used  to  render  [Af]  and  [K]  diagonal  is  a  transformation  to 
normal  coordinates.  The  problem  is  to  find  an  orthogonal  matrix*  d>  which 
simultaneously  diagonalizes  [Af]  and  [AT].  If  [<I>]  is  orthogonal  then  for  unity  modal 
mass  normalization 

[<S>f[M]m^[I]  ,  (3.43) 


or 


[<i>]^[Af  ]  = 

Recall  that  the  identity  matrix  can  be  defined  as 

0*^4^  *  0^0  “  1 


where 

10.. 
0  10. 
^  .  0  1  0 
..01 


(3.44) 


(3.45) 


The  process  of  simultaneously  diagonalizing  the  [Af]  and  [K]  matrices  can 
begin  by  rearranging  the  eigenproblem  equation  (3.35)  as 

[An$  =  a)*(Af]$  .  (3.46) 

Corresponding  to  the  ith  and  jth  normal  modes  are  the  vector  equations 

[K]^^  *  iohMh'  .  (3.47) 

and 

[K]^  =  .  (3.48) 

Using  the  identity  matrix  to  make  them  scalar  equations,  equation  (3.47)  is 

•T 

multiplied  by  ^  ,  where  the  transpose  changes  the  n  x  1  column  matrix  into  a 
1  X  n  row  matrix,  and  equation  (3.48)  is  multiplied  by  so  that 

=(OiV’'[Af]<^‘  ,  (3.49) 
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and 


.  (3.50) 

Since  [K]  and  [M]  are  real,  symmetric  matrices  [X"]  =  [A’]^  and  [Af]  =  [M]^. 
Making  use  of  this  property,  the  transpose  of  equation  (3.50)  yields 

.  (3.51) 

which  when  subtracted  from  equation  (3.49)  gives  the  scalar  equation 

0  =  ((i)i*-a)/)^’‘lM)$‘  .  (3.52) 


If  (Oi  ,  there  are  two  different  normal  modes,  then 


♦i  W 


^11  ^21 

Af  21  Af  22 


♦  l 

♦I 


=  0 


(3.53) 


One  of  the  key  characteristics  of  normal  modes  is  that  every  normal  mode  is 
orthogonal  to  every  other  mode.  Specifically, 


r 


i  *j 
i  • 


(3.54) 


The  mass  matrix  [M]  is  now  diagonalized  into  a  new  modal  mass  matrix  M^.  This 
matrix  M}  is  a  measure  of  how  much  mass  is  participating  in  the  ith  normal  mode. 
Similarly  for  the  non-diagonal  stiffiiess  matrix  [K] 


= 


Jo  i 

K,  i  =J 


(3.55) 


Therefore,  diagonalization  is  achieved  by  a  linear  transformation  using  the 
orthogonal  modal  matrix  <1>,  now  defined  as 


[<!>]  = 


W]  (♦'  Jf 


(3.56) 


which  represents  the  spatial  and  temporal  parts  of  the  mode  shape.  It  is  composed 
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of  the  eigenvectors  and  V.  The  eigenvalues  form  a  spectral  matrix 


0 

0  (0* 


Returning  to  the  equation  of  motion 

[Af]x  =0  .  (3.35) 

and  introducing  the  orthogonal  modal  matrix  as  part  of  the  solution 

X  =  [<!>]£  .  (3.57) 

which  when  substituted  back  into  equation  (3.35)  yields 

[MKOlj =0  .  (3.58) 

To  perform  a  linear  transformation,  this  equation  is  multiplied  by  [<I>]^  so  that 

=0  .  (3.59) 

which  diagonalizes  [M]  and  [F  \  to  create  a  new  normalized  relation 


Mu  0  .  .  ' 

I 

0 

0  M22  0 

r  • 

0  K22  0 

0  .  0 

0  M(j 

[9 

+ 

0  .  0 

0  Ky 

=  0 


(3.60) 


Now  the  n  coupled  differential  equations  are  uncoupled.  Every  ith  normal  mode 
now  has  its  own  ith  row  which  can  be  lifted  out,  such  as 


M|7i  +  Kj(7i  =  0 

which  is  like  the  familiar  harmonic  oscillator  mass/spring  system  where 

•  ^  n 


(3.61) 


or 


Qi  +  ^<li  =  0 
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Any  arbitrary  motion  of  the  system  can  thus  be  written  as  a  linear  combination  of 
the  eigenvectors 

=  9i[«i] +<72[«2]  .  (3.62) 

where  Qi-  Acos (Oit  and  <12-  Acos .  Here  A  is  an  amplitude  coefEcient  based 
on  the  perturbing  force.  Hence,  the  finite  element  computer  code  solves  for  the 
normal  modes  of  a  structure  by  finding  a  modal  matrix  <I>  which  will  diagonalize 
[ilf  ]  and  [Kl 

Therefore  a  normal  mode  is  made  up  of  a  particular  eigenvalue  (frequency) 
and  a  particular  eigenvector  (mode  shape),  which  is  orthogonal  to  every  other 
normal  mode.  A  normal  mode  describes  how  much  mass  and  stif&iess  is  involved 
in  the  structure  as  it  vibrates.  That  is,  the  coordinate  is  a  measure  of  how  much 
the  mode  shape  d>,  is  present  The  value  of  finding  the  normal  modes  of  a 
structure  is  that  any  oomplez  motion  can  be  represented  as  a  linear  combination  of 
normal  modes.  For  such  an  intricate  system  as  the  PHALANX  gun,  where  there 
are  many  components  contributing  to  its  dynamic  response,  to  find  its  normal 
modes  suddenly  makes  it  a  far  easier  system  to  understand.  That  is,  if  the  normal 
modes  of  the  gun  are  obtained,  which  then  forms  a  modal  component  model,  when 
excitation  forces  are  applied,  the  resulting  displacements  are  the  result  of  exciting 
certain  normal  modes.  An  analogy  to  this  is  that  any  color  can  be  made  up  of 
various  amounts  of  the  primary  colors  red,  blue  and  yellow.  So  in  order  to  obtain  a 
set  of  accurate  normal  modes  the  finite  element  model  must  be  as  close  a 
representation  of  the  real  gun  as  possible  which  is  the  topic  of  the  next  chapter. 


•*1 

Jt2 
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IV.  MODEL  DESCRIPnON 


A.  TYPES  OF  FINITE  ELEMENTS  USED  IN  THE  MODEL 

At  the  outset  of  this  project,  the  gun  barrels  alone  were  thought  to  be  a 
primary  source  of  dispersion.  In  fact,  previous  to  this  project,  the  Navy  had 
already  let  a  contract  for  a  new  barrel  design.  So,  the  initial  proposal  for  research 
by  the  Naval  Postgraduate  School  called  for  a  finite  element  model  of  the  six-barrel 
system.  However,  in  consultation  with  Mike  Hatch,  a  private-industry  finite 
element  expert,  it  was  agreed  that  in  order  to  truly  represent  the  complex 
dynamics  of  the  gim,  it  would  be  important  to  model  not  just  the  barrels,  but  the 
whole  gun  assembly  through  to  its  grounding  points  to  the  g\in  mount.  This 
holistic  approach,  while  greatly  expanding  the  scope  of  the  work,  would  not  only 
offer  a  much  better  assessment  of  the  g\in  system,  but  would  greatly  expand  its 
applications  as  well.  The  school  would  also  be  the  first  organization  to  model  the 
complete  gun  assembly,  thxis  offering  a  complete  model  in  one  package.  The 
benefits  of  such  a  model  are  many.  Once  the  model  was  assembled  and  the 
dynamic  response  validated  against  the  current  dispersion  pattern,  new 
components  could  be  "swapped  out"  in  the  model  and  their  effect  on  dispersion 
would  provide  a  valuable  assessment.  Not  only  could  individual  components  such 
as  the  barrels  be  tested,  but  combinations  of  components  could  be  tested  together. 
This  model  therefore,  would  go  a  long  way  toward  understanding  the  dynamics  of 
the  gun  as  well  as  providing  a  valuable  tool  for  finding  ways  to  tighten  the 
dispersion  pattern. 

Once  blueprints  of  the  model  were  received,  the  next  big  decision  involved  the 
appropriate  selection  of  finite  elements  to  represent  each  component.  In  all  there 
are  ten  different  parts  of  the  gun,  coimting  duplications  such  as  the  barrels,  as  just 
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one  component.  The  finite  element  modeling  experience  Professor  Gordis  and 
Mike  Hatch  was  helpful  in  deciding  which  areas  of  the  gun  system  should  be 
modeled  with  great  detail  and  areas  where  time-consuming  attention  could  be 
safely  passed  over.  Table  4-1  lists  the  components  and  the  corresponding  element 
choices.  Before  discussing  the  properties  of  the  elements  that  were  used  in  the 
model,  it  is  important  to  give  an  overaU  view  on  element  selection.  Each  element 
comes  with  a  price.  Some,  such  as  bricks,  are  more  expensive  to  use  than  others 
since  they  require  more  nodal  points  to  define  them,  in  this  case  eight.  Others, 
such  as  beam  elements,  as  discussed  in  the  previous  chapter,  require  only  two 
nodes.  In  order  to  shorten  the  computational  time  and  make  the  resulting  modal 
shape  animation  as  smooth  as  possible,  it  is  important  to  choose  as  many  less 
expensive  elements  as  possible.  At  the  same  time,  to  accurately  represent  the  real 
gun,  the  expensive  elements  are  sometimes  the  best  choice.  Hence,  it  is  a  balance 
that  often  goes  through  several  iterations  before  a  final  choice  is  made.  As  can  be 
seen  in  Table  4-1,  the  gun  is  comprised  of  five  categories  of  elements,  with  many 
different  variations  within  each  category,  such  as  changes  in  dimensions,  material 
properties,  etc.  The  following  subsections  describe  the  general  properties  of  each 
category  of  element  used  in  the  model. 

1.  Brick  Elements 

As  mentioned  before,  the  solid  brick  elements  are  the  most  complex 
elements  to  use  in  that  they  require  eight  nodes  to  define.  A  single  brick  element 
is  shown  in  Figure  4-1.  A  brick  element  is  one  of  the  most  versatile  elements  to 
use  since  it  allows  modeling  of  any  structure  with  any  geometry  [Ref  10].  The 
general  geometry  of  the  brick  is  that  its  width  is  roughly  comparable  to  its  length. 
It  is  a  three-dimensional  solid  whose  nodes  are  allowed  to  translate  along  each 
axis,  but  not  rotate.  Therefore,  it  has  only  three  degrees  of  fireedom.  It  is  best 
used  to  represent  solid,  invariant  structures  without  many  intricate  details  or 
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curvature,  such  as  engine  blocks  or  in  our  case  the  breech.  Fvurthermore,  more 
bricks  are  required  to  reflect  bending  and  shear  than  do  thin  shells  or  beam 
elements.  A  structure  represented  with  bricks  often  looks  the  most  realistic. 
However,  it  comes  with  the  price  of  reduced  accuracy,  especially  if  not  enough 
bricks  are  used,  and  in  any  case,  it  is  usually  the  most  inefficient  element  to  use 
[Ref.  11]. 


TABLE  4-1.  ELEMENT  SELECTION 


Elementary  Model 

r— - 

Component 

Element  Types 

Muzzle  Clamps 

Bricks,  Thin  Shells,  Rigid  Elements 

Mid-Barrel  Clamps 

Bricks,  Rigid  Elements 

Barrels 

Beams  with  tubiilar  cross-section 

Stub  Rotor 

Bricks,  Thin  Shells 

Rotor 

Bricks 

Needle  Bearing 

Bricks,  Springs 

Casing 

Thin  Shells 

Thrust  Bearing 

Bricks,  Springs 

Recoil  Adapters 

Bricks,  Rigid  Elements,  Springs 

Ball  Joint 

Rigid  Elements,  Springs 
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Figure  4*1.  A  solid  linear  brick  element. 


2.  Thin  Shell  Elements 

A  thin  shell  or  plate  element  is  a  two^limensional  solid  whose  thickness 
is  small  compared  to  its  length  or  in-plane  dimensions.  A  thin  shell  element  is 
shown  in  Figure  4-2.  It  is  a  much  more  efficient  element  to  use  than  a  solid  brick 
since  it  only  requires  four  nodes  to  define  its  shape.  Its  thickness  is  defined  when 
the  element  is  created,  but  it  does  not  show  up  in  the  display.  The  best 
applications  for  it  are  in  modeling  situations  where  a  wall  or  membrane  has  a 
nearly  constant  thickness  or  its  thickness  varies  linearly  from  one  edge  to  another. 
The  element’s  nodes  are  allowed  to  not  only  translate,  but  rotate  as  well,  thus 
giving  it  six  degrees  of  fireedom.  These  additional  degrees  of  fireedom  allow  fewer 
thin  shell  elements  to  be  used  than  brick  elements  in  obtaining  the  same  accuracy 
for  bending  and  shear  motion.  One  disadvantage  in  using  thin  shell  elements  is 
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that  they  are  not  able  to  model  stresaes  that  vary  through  ita  thickness  as  can 
brick  elements  [Ref.  11]. 


Figure  4*2.  A  thin  shell  element. 


3.  Beam  Elements 

Beam  elements  are  the  simplest  elements  in  that  they  only  require  two 
nodes  to  define  them.  The  cross  section  of  the  beam  is  defined  when  the  beam  is 
created  and  can  have  as  detailed  and  complex  a  shape  as  the  user  desires.  Figure 
4*3  shows  a  beam  element  as  it  appears  in  the  display  of  the  finite  element  model 
and  its  cross  section  which  only  appears  when  referencing  the  beam’s  properties. 
This  cross  section  is  invariant  along  the  element’s  length.  To  represent  a  beam  of 
varying  cross-section,  several  elements  are  reqxiired.  For  example,  the  gun  barrels 
are  represented  using  beam  elements,  where  the  barrels’  varying  thickness  is 
reflected  in  the  varying  cross-section  of  12  elements.  The  nodes  are  allowed  the 
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full  six  degrees  of  fireedom  which  is  usehil  in  displaying  bending  and  deflection. 
Again  referring  to  the  modeled  gun  barrels,  the  beam  elements  provide  a  smooth 
animation  picture  which  gives  better  insight  into  the  dynamics  of  the  gun.  A 
possible  problem  with  using  beam  elements  is  that  they  do  not  show  stress  through 
their  cross  section  and  the  cross  sections  can  be  easily  oriented  the  wrong  way  [Ref. 
11].  However,  in  most  cases,  there  is  more  interest  in  the  maximum,  outer  fiber 
stress  than  that  of  the  inner  stress. 


Figure  4-3.  (a)  A  beam  element,  (b)  A  representative  beam  element’s 
cross  section. 


4.  Spring  Elements 

Spring  elements  pan  be  defined  with  two  nodes  (node*to*node)  or  in  terms 
of  one  node  (node-to-ground).  They  can  also  be  either  a  translational  spring 
operating  along  a  defined  axis,  or  a  rotational  spring  which  represents  the  bending 
moment  or  torque  that  is  required  to  rotate  one  node  relative  to  another  or  to 
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ground  [Ref.  12].  Figure  4-4  shows  a  node-to-node  and  a  node-to-ground  spring. 
The  stifhiess  is  defined  along  each  axis  or  un.^udally  for  a  translational  spring  or 
about  an  axis  for  a  rotating  spring.  Three  rotational  degrees  of  fireedom  are 
assigned  to  each  node.  In  the  PHALANX  gun  model  they  were  used  to  represent 
the  stiflhess  of  the  ball  bearings  and  the  recoil  shock  absorbers. 


m - * 


Figure  4-4.  (a)  A  node-to-node  spring  element,  (b)  A  node-to-ground 
spring  element. 


5.  Rigid  Elements 

Rigid  elements  are  massless,  infinitely  stiif  bars  used  to  connect 
components  together.  That  is,  the  rigid  element  implies  that  the  motion  of  the 
nodes  on  the  element  are  to  be  related  to  each  other  as  if  they  were  connected 
together  by  an  infinitely  rigid  bar  [Ref.  12].  Two  nodes  that  are  connected 
together  with  a  rigid  element  move  together  according  to  a  linear  algebraic 
expression  called  a  constraint  relation.  An  example  of  such  a  relation  is  given 
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below. 


11-0.5x2  =  0  (4.1) 

In  this  example,  the  first  term  describes  the  distance  node  1  moves  in  the  x 
direction  and  the  second  term  that  of  node  2  in  the  z  direction.  The  right  hand 
side  of  the  equation,  which  is  a  constant  term,  must  be  zero  when  solving  for  the 
normal  modes  of  vibration.  This  equation  can  be  written  in  matrix  notation  as 


0.5 

*2 

.^2. 

1 

.*2 

Representing  the  coefficient  matrix  with  T,  this  relation  can  then  be  written  as 

^original  =  »  (4.3) 

where  the  nodal  position  is  defined  in  terms  of  a  coefficient  matrix  and  the 
constrained  node  matrix  Xg .  Substituting  thi*  equation  into  the  equation  of  motion 

F^KTx^  ,  (4.4) 

and  multiplying  by  the  transpose  of  the  coefihcient  matrix  T  yields 

TJF  =  r/XTxc  ,  (4.5) 

which  reduces  to  the  constrained  equation  of  motion, 

Fc=K^x^  .  (4.6) 

This  is  an  example  of  how  constraint  equations  transform  an  unconstrained 
stiffiiess  matrix  K  to  a  constrained  stiffness  matrix  Kg.  A  rigid  element  can  be 
\ised  to  connect  firom  two  to  thirty-two  nodes  together.  Figure  4-5  shows  a  rigid 
element  connecting  five  nodes  together.  Each  node  at  the  element’s  end  points  are 
allowed  all  translations  and  rotations  that  are  allowed  for  that  node.  For  example, 
if  the  rigid  element  is  connected  to  a  node  of  a  brick  element,  then  it  is  only 
allowed  to  translate.  In  the  PHALANX  gun  model,  rigid  elements  were  used  to 
connect  the  barrels  to  the  clamps  and  stub  rotor  and  the  recoil  adapter  springs  to 
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the  g\in  body.  "End  release  codes"  can  also  be  defined  allowing  the  simulation  of 
hinges  or  ball-and-socket  joints  [Ref.  12]. 


Figure  4-5.  A  rigid  element  connecting  five  nodes  together. 

A  further  classification  within  each  element  category  is  the  element  order. 
That  is,  the  number  of  nodes  per  element  can  be  increased  to  increase  the  order  of 
the  equations  used  to  interpolate  the  strain  between  them.  There  are  two  nodes 
required  for  each  edge  in  a  linear  element,  three  for  a  parabolic  and  four  nodes  for 
a  cubic  [Ref.  11].  All  the  elements  used  in  this  model  have  linear  order. 

B.  CONSTRUCTION  OF  THE  PHALANX  GUN  MODEL 

A  breakdown  of  the  components  used  in  the  construction  of  the  PHALANX 
gun  model  is  shown  in  Figure  4-6.  This  is  a  hidden  line  image  which  shows  the 
geometry  of  the  elements  on  the  side  of  the  gun  parts  facing  the  viewer.  Figure  4-7 
is  a  shaded  image  of  the  completely  assembled  gun.  Here,  the  individual  elements 
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cannot  be  seen,  but  the  effect  of  lighting  and  solid  colors  gives  a  realistic  picture  of 
the  gun  it  represents.  Both  figures  provide  a  brief  overview  of  how  the  gun  is 
assembled.  The  six  barrels  each  slide  through  the  stub  rotor  and  are  screwed  into 
the  rotor  a  distance  of  approximately  7.6  cm  (3  inches).  The  stub  rotor  itself  is 
bolted  to  the  rotor  and  serves  to  support  the  barrels  near  their  base.  The  six- 
barrels  are  connected  together  with  a  mid-barrel  clamp  and  a  forward  muzzle 
clamp.  This  whole  assembly  is  now  attached  to  the  gun  body  through  a  rear  needle 
bearing  and  a  double-roller  contact  bearing.  The  assembled  gxm  body  is  coupled  to 
the  rest  of  the  CIWS  mount  at  three  locations.  Two  recoil  adapters  are  positioned 
on  either  side  of  the  gun  body  where  they  connect  the  gun  body  to  the  mount  via  a 
pivot  joint.  The  third  grounding  point  is  a  ball  joint  at  the  rear  of  the  gun  body. 
All  of  these  components  will  be  explained  in  further  detail  in  the  following 
subsections. 

1.  Gun  Barrels 

The  six  rifled  barrels  are  approximately  150  cm  (59  inches)  long  with  an 
inner  diameter  of  2  cm  (0.787  inches)  and  weigh  about  8.2  kg  (18  lb).  The  barrels 
are  secured  to  the  breech  with  interrupted  locking  lugs  that  mate  after  the  barrels 
are  screwed  in  past  three  threads  [Ref.  13}.  As  mentioned  before,  the  six  barrels 
are  represented  using  beam  elements.  Each  barrel  is  made  up  of  39  elements  with 
12  different  cross-sections.  Figure  4-8  shows  a  single  barrel  with  its  elements 
separated. 

2.  Muzzle  Clamp 

The  muzzle  clamp  restricts  barrel  movement  near  the  barrel  tips.  It 
consists  of  two  plates  held  together  with  a  steel  web  meshing.  The  clamp  also 
restrains  individual  barrel  movement  during  firing  and  a  self-locking  nut  is  used  to 
secure  the  clamp  to  all  the  barrels  [Ref.  13].  The  plates  are  represented  using 
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Mac  support 


Figure  4-6.  A  hidden  line  image  of  the  diaassembled  components  of  the 
PHALANX  gun. 

brick  elements  while  the  web  meshing  is  modeled  with  thin  shells  as  can  be  seen  in 
Figure  4-9. 

3.  Mid-Barrel  Clamp 

The  mid-barrel  clamp  is  positioned  roughly  near  the  center  of  the  barrel 
group.  It  physically  locks  the  barrels  to  the  clamp  with  tabs  that  fit  into  cutouts  in 
each  barrel.  It  prevents  the  barrels  from  becoming  unscrewed  in  the  rotor  during 
firing.  The  mid-barrel  clamp,  as  shown  in  Figure  4-10,  is  represented  entirely  with 
bricks. 
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Figure  4*7.  A  shaded  image  of  the  PHALANX  gun. 


4.  Stub  Rotor 

The  stub  rotor  which  is  bolted  to  the  rotor  supports  the  base  end  of  the 
six  barrels.  It  permits  alignment  of  the  barrels  so  that  the  locking  lugs  ran  be 
engaged  [Ref.  13].  like  the  muzzle  clamp  it  consists  of  two  plates  separated  by  a 
metal  meshing.  Figure  4-11  shows  the  modeling  of  it  with  bricks  and  thin  shells. 

5.  Rotor 

The  rotor  is  the  most  massive  component  in  the  gun  body.  It  is 
hydraulically  rotated  at  a  speed  of  750  rpm  (12.5  Hz)  in  the  coimterclockwise 
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Figure  4-8.  A  single  barrel  with  with  39  beam  elements. 

direction  when  viewed  firom  the  rear  [Ref.  13].  The  gun  assembly,  which  is  shown 
in  Figure  4-12,  is  bolted  onto  the  front  of  it.  The  rotor  is  shown  in  FHgiire  4-13  and 
is  made  entirely  of  bricks.  The  rotor  is  coupled  to  the  gun  body  through  the  rear 
needle  bearing  and  the  forward  double  roller  contact  bearing. 


Figure  4*9.  The  muzzle  clamp  represented  with  brick  and  thin  shell 
elements. 

6.  Double  Row  Contact  Bearing 

The  double  row  angular  contact  bearing,  often  referred  to  in  this  text  as  a 
thrust  bearing,  consists  of  two  rows  of  ball  bearings  for  a  total  of  76  balls.  The 
inner  race  attaches  to  the  rotor  while  its  outer  race  attaches  to  the  gun  body.  Its 
purpose  is  to  align  the  rotor/gun  assembly  and  transfer  the  firing  load  to  the  rest  of 
the  g\in  body.  It  provides  both  radial  and  axial  support.  Figure  4-14  shows  the 
races  represented  by  bricks.  The  ball  bearings  are  modeled  using  springs  with  the 
loaded  contact  angle  being  represented  by  an  angle  determined  through  an 
accepted  computer  program  from  Mike  Hatch. 
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Figure  4-10.  The  mid-barrel  clamp  represented  with  bricks. 

After  the  decision  was  made  to  model  the  whole  gun  subsystem,  careful 
attention  was  given  to  the  blueprints  of  the  gun  body.  It  then  became  readily 
apparent  that  the  two  bearings  were  playing  a  critical  role  in  the  dj'namics  of  the 
gun.  They  provide  the  only  means  of  support  to  the  rotating  and  firing  six-barrel 
gun  assembly.  Therefore,  the  intricate  details  of  the  bearings,  especially  the  double 
row  bearing,  were  modeled  as  can  be  seen  firom  a  cutaway  view  showing  the  cross 
section  of  the  double  row  bearing  in  Figure  4-15.  It  is  also  noteworthy  that  this 
bearing  is  mounted  to  the  rotor  through  a  narrow  lip  which  was  also  accurately 
modeled  and  is  pointed  out  in  that  figure.  The  entire  firing  and  rotating  forces  are 
transmitted  through  this  lip  which  seems  to  be  subjecting  it  to  a  severe  stress. 
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Figure  4-11.  The  stub  rotor  represented  with  bricks  and  thin  shells. 

Obtaining  the  axial,  radial  and  torsional  stiilhess  values  for  the  double 
row  contact  bearing  was  difficult  since  the  manufacturer,  Kaydon  Bearing  in 
Muskegon,  Michigan  did  not  have  records  for  this  bearing  in  this  application. 
However,  after  supplying  Kurt  Sheridan  of  Kaydon  with  the  reaction  force 
information,  determined  from  the  finite  element  model,  he  was  able  to  provide 
approximate  figures  of  the  bearing  stiffness  and  contact  angle.  With  this 
information  Mike  Hatch  was  then  able  to  provide  spring  stiffness  values  and  an 
angle  with  which  to  mount  the  springs  in  relation  to  the  inner  and  outer  races. 
There  are  a  total  of  36  springs,  one  for  each  of  the  18  elements  used  to  define  its 
circumference  then  multiplied  by  two  to  represent  the  double  row  of  bearings. 
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Figure  4-12.  The  gun  barrels  with  clamps  and  stub  rotor. 

7.  Needle  Bearing 

The  needle  bearing  provides  the  rear  support  for  the  rotor.  It  serves  only 
as  a  radial  support  and  not  as  a  thrust  absorber.  The  inner  race  of  the  bearing  is 
attached  to  the  rear  wall  of  the  gun  body  while  the  outer  race  is  attached  to  the 
rotor.  Figure  4-16  is  a  translucent  picture  showing  the  inner  assembly  of  the  gun 
body.  In  this  figure  it  is  possible  to  see  both  the  needle  bearing  and  the  double  row 
bearing  and  their  relation  to  the  rotor  and  gun  body.  Similar  to  the  double  row 
bearing,  the  bearing  stiffiiess  is  represented  with  18  springs.  In  order  to  correctly 
represent  the  actual  bearing  stifihess,  ^6(aniv>  effective  spring  stiffness,  K^ff, 
must  be  calculated.  No  matter  which  way  the  rotor  moves  on  the  real  gun,  it 
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Figure  4-13.  The  rotor  represented  with  bricks. 


experiences  a  constant  bearing  stif&iess,  whereas  in  the  modeled  gun  which  has 
only  18  springs,  the  angle  between  them  must  be  taken  into  account.  Figure  4-17 
demonstrates  the  geometry  needed  to  make  this  calculation.  Consider  a  force 
applied  in  any  one  direction  to  the  model  bearing.  The  equation  of  motion  would 
be 


F^Ku^X 


(4.7) 


Here,  Ktot  is  the  total  spring  stifEaess  value  that  would  take  into  accoiint  all  18 
springs.  However,  Ktu  *  F^taring  since  as  it  stands  now,  the  modeled  bearing  will 
present  a  larger  spring  stiffness  than  the  actual  bearing.  To  correct  for  this 
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Figure  4-14.  The  double  row  contact  bearing. 


problem,  the  20P  angle  between  the  springy  must  be  described  by  the  relation 

Ktot  =  i,K,ffCoaai  .  (4.8) 

>>o 


Solving  for  Kgff  gives 

t-18 

JCgf  f  ^  COSOtj  —  '^bcarinf  ' 
i-0 

where 


^bearing  =Kgffi2c08(f  +  4cO8  20®  +  4c08A(f  +  .  .  .) 


Rearranging  terms  gives, 


(4.9) 


(4.10) 
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Figure  4>15.  Cross  sectional  view  of  the  double  row  contact  bearing. 
Note  the  narrow  lip  which  attaches  to  the  rotor. 


= 


^bearing 

11.5175 


(4.11) 


Therefore,  after  obtaining  the  stifhaess  value  of  the  needle  bearing  from  the 
manufacturer,  it  was  divided  by  a  factor  of  11.5175  to  fit  the  model  [Ref.  14]. 


8.  Gun  Body 

The  gun  body  is  a  steel  housing  which  contains  the  bearings,  rotor  and 
stub  rotor  as  shown  in  Figiure  4*16.  Since  the  wall  thickness  is  only  0.5  cm  (0.2 
inches)  and  is  small  compared  to  its  other  dimensions,  it  is  represented  by  thin 
shells.  It  is  part  of  the  load  path  between  the  firing  barrels  and  ground.  On  either 
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Figure  4-16.  Translucent  image  of  the  gun  body.  It  shows  the  relation 

between  the  two  bearings  and  the  rotor.  Also  visible  are  the 
recoil  adapters  and  the  rear  support. 

side  of  the  gun  body,  roughly  opposite  to  the  double  row  bearing,  are  the  two  recoil 
adapters.  At  the  rear  of  the  gun  body  is  the  third  attachment  to  ground,  the  ball 
joint. 

9.  Recoil  Adapters 

The  two  recoil  adapters  serve  to  absorb  a  portion  of  the  forces  caused 
when  the  gun  is  fired  and  then  the  associated  counter-recoil  force  [Ref  13].  Figure 
4-18  is  a  photograph  of  the  PHALANX  gun  on  the  USS  ARKANSAS  which  shows 
the  construction  and  mounting  of  an  actual  recoil  adapter.  Since  the  overall  shape 
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Figure  4*17.  Spring  placement  for  the  needle  bearing. 

of  the  recoil  adapter  plays  a  minor  role  in  the  dynamics  of  the  gun,  it  was  possible 
to  represent  its  recoil  spring  stifhiess,  placement  and  attachment  to  ground  simply, 
but  accurately.  Looking  at  Figure  4*16,  the  recoil  adapter’s  mass  and  stiffness  was 
modeled  using  brick  elements  in  the  form  of  "ears”  sticking  out  on  either  side  of 
the  gun  body.  In  order  to  get  the  correct  placement  of  the  roller,  rigid  elements 
were  xised  to  coimect  the  brick  element  ears  to  the  recoil  spring.  In  turn,  the  recoil 
spring  was  connected  to  a  roller.  The  specifics  of  the  roller  will  be  discussed  in  the 
bovmdary  conditions  section.  The  stiffness  value  for  the  spring  was  obtained  from 
a  hysteresis  curve  as  shown  in  Figure  4-19  [Ref.  15].  Here  the  slope  for  the  upper 
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and  lower  branches  of  the  curves  was  calculated  and  averaged  to  get  the  stiflhess 
value  \ised  in  the  model. 


Figure  4-18.  Photograph  of  the  gun  body.  Note  the  arrangement  of  the 
recoil  adapters  and  their  construction. 


10.  Ball  Joint 

The  rear  ball  joint  provides  the  third  support  of  the  tripod  supported  gun. 
The  geometry  of  it  is  displayed  in  Figure  4*20.  It  consists  of  a  3.84  cm  (1.5  inch) 
diameter  ball  supported  by  a  9.63  cm  (3.79  inches)  long  rod  2.54  cm  (1  inch)  in 
diameter  which  in  turn  is  attached  to  the  gun  mount.  The  ball  fits  into  a  socket 
mounted  on  the  back  of  the  rear  wall  of  the  gun  body.  When  the  gun  fires,  it  sets 
back  by  about  half  an  inch  and  oscillates  about  the  new  equilibrium  point.  The  ball 
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Figure  4-19.  Hysteresis  cxirve  of  the  recoil  adapter  stiflhess. 

slides  into  the  socket  which  only  provides  radial  support.  Since  this  is  a  relatively 
simple  arrangement,  a  separate  side  model  was  built  to  establish  the  stiilhess  of 
the  socket.  This  stiffiiess  was  represented  by  two  springs  along  the  y  and  z  axes 
perpendicular  to  the  length  of  the  gun  thus  providing  the  radial  support.  The 
springs  were  attached  to  groimd  on  one  end  while  their  other  end  connected  to 
rigid  elements  that  attached  to  the  gun  body.  These  rigid  elements  represented 
the  support  given  by  the  cup  to  the  gun  body.  Figure  4-21  shows  this  arrangement 
in  the  computer  model. 
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Figure  4*20.  Blueprint  of  the  rear  ball  joint. 

C.  BOUNDARY  CONDITIONS 

Having  modeled  the  components  of  the  gun,  the  next  msyor  step  is  to  ensure 
that  the  proper  boimdary  conditions  are  reflected  in  the  model.  The  boundary 
conditions  define  how  the  components  are  connected  to  each  other  as  well  as  how 
the  whole  g\m  is  attached  to  the  gun  mount.  The  next  category  of  boundary 
conditions  that  must  be  determined  are  the  types  of  degrees  of  freedom. 

As  the  model  stands  now,  it  represents  an  idealized  gun  with  no  worn 
bearings,  or  recoil  springs  and  no  tolerance  between  connected  components.  It  may 
be  better  than  a  brand  new  gun.  As  mentioned  earlier,  the  barrels  are  attached  to 
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Fissure  4*21.  Representation  of  the  ball  joint  in  the  model. 

the  clamps,  stub  rotor  and  rotor  with  rigid  elements.  Figure  4*22  shows  an 
example  of  a  barrel  connected  to  a  clamp.  The  stub  rotor  is  connected  to  the  rotor 
with  rigid  elements  as  well  The  rotor  is  connected  to  the  gun  body  through  the 
two  bearings.  The  gun  body  in  turn  is  connected  to  ground  via  the  recoil  adapters 
and  ball  joint. 

Figure  4*16  expands  further  on  the  specific  ranges  of  motion  allowed  at  these 
three  supporting  points  showing  how  they  were  defined  in  the  model.  At  the  two 
pins  which  hold  the  recoil  adapters  to  ground  in  the  real  gun,  they  are  allowed  to 
rotate  about  the  z  or  horizontal  axis,  but  do  not  rotate  or  translate  in  any  other 
direction.  To  properly  represent  this  condition  in  the  real  gun,  the  node  between 
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Figure  4*22.  Example  of  rigid  elements  connecting  a  barrel  to  a  clamp. 

the  rigid  elements  and  the  springs  is  allowed  to  move  in  the  z  or  axial  direction 
and  to  rotate  about  the  z  axis.  The  translation  in  the  x  direction  is  necessary  in 
order  to  allow  the  springs  to  be  exercised.  The  springs  in  turn  are  attached  to 
groimd  in  order  to  satisfy  the  real  boundary  conditions,  albeit  broken  up  into  two 
locations,  given  the  unique  representation  of  the  recoil  adapters  in  the  model. 
Motion  at  the  ball  joint  of  the  real  gun  is  in  the  fore-all  or  x  direction,  with  the  cup 
providing  restraint  for  translational  motion  in  the  y  or  z  directions.  All  rotations 
are  allowed,  however  due  to  the  two  forward  supports  this  motion  would  be 
unlikely.  In  the  computer  model,  the  springs  restrict  motion  in  the  y  and  z 
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directions,  thus  like  the  cup,  and  there  are  no  restrictions  for  motion  in  the  x 
direction. 

Unlike  the  real  gun,  this  model  does  not  spin.  This  is  where  kinematic  and 
unknown  degrees  of  freedom  come  into  play.  The  normal  mode  of  rotation  about 
the  X  axis  is  prevented.  Since  the  normal  modes  are  being  solved  without 
gyroscopics,  that  is,  without  an  unbalance  (centrifugal)  excitation,  it  is  restricted  in 
the  model  by  using  kinematic  degrees  of  freedom.  Four  kinematic  degrees  of 
freedom  were  used  to  restrict  the  rotation  of  the  rotor  inside  the  gun  body.  This 
does  not  substantially  degrade  the  understanding  of  the  gun’s  dynamic  motion, 
since  the  excitation  force  occurs  at  the  same  place  and  is  virtually  independent  of 
the  rotation.  That  is,  as  the  rotor  is  turning  and  the  barrels  are  sequentially  firing, 
they  all  fire  at  the  same  location.  The  barrels  are  not  actually  rotating  in  the 
computer  model,  but  the  excitation  force  is  applied  at  the  base  of  a  barrel  which  is 
in  the  proper  location  to  represent  each  barrel  as  it  fires  in  the  rotating  gun.  The 
unknown  degrees  of  freedom  correspond  to  the  differential  equations  solved  by  the 
computer.  In  this  model,  there  are  22,470  degrees  of  freedom. 

D.  COMPARISON  BETWEEN  THE  REAL  AND  MODELED  GUN 

After  completion  of  the  PHALANX  gun  model,  one  of  the  first  steps  toward 
validation  is  to  compare  its  properties  with  those  of  the  real  gun.  Based  on 
information  provided  by  Dick  Sirola,  the  properties  used  for  comparison  were  mass, 
center  of  gravity,  and  moments  of  inertia  for  each  of  the  three  coordinate  axes. 
This  information  and  their  percent  difference  between  the  computer  model  and  the 
real  gvm  is  presented  in  Table  4'2.  Fturther  comparisons  can  be  made  and  the 
model  refined  to  provide  a  better  match.  These  small  differences  are  not  important 
to  the  conclusions  reached  by  this  study. 
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TABLE  4-2.  COMPARISON  BETWEEN  REAL  AND  MODELED  GUN 


Comparison 


Mua  (kg)  Center  of  Gravity  cm)  Moment  of  Inertia  (kg-m 


EM 

OH 

/» 

tyy 

— 

M61A1 

113.5 

53.83 

0.04 

0.0 

6796 

2.38e® 

2.38e® 

Model 

120.5 

55.2 

-0.304 

-0.053 

7240 

2.59e® 

HHIHH 

%Diff. 

5.81 

2.48 

5.75 

5.40 

6.13 

8.01 

8.40 

Another  important  area  of  oomparison  is  to  determine  whether  the  actual 
forces  in  the  gun  are  linear  or  non-linear.  The  computer  will  only  solve  the  linear 
problem.  Experimental  data  verifies  that  the  stifihess  calculations  are  linear  for 
the  range  of  barrel  tip  motion  less  than  -  0.25  cm  (0.1  inch)  [Ref.  16].  Figure  4*23 
shows  the  test  results  of  muzzle  deflection  versus  a  load  up  to  300  lb  applied  in  the 
vertical  direction  at  the  barrel  tip.  Since  the  dispersion  data  shows  that  barrel  tip 
displacement  is  less  than  =  0.25  cm,  the  system  is  considered  linear.  The  linear 
computer  solution  is  compatible  with  the  experimental  data. 
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Figure  4-23.  Experimental  data  of  muzzle  deflection  versus  load. 
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V.  MODES  OF  VIBRATION 


A.  MODES  OF  VIBRATION  OF  THE  GUN  BARRELS 

One  of  the  first  components  of  the  model  to  be  developed  were  the 
barrels.  Two  different  levels  of  modeling  detail  were  used  to  create  two  different 
versions  of  the  barrel.  The  first,  called  the  advanced  barrel,  involved  the  use  of 
496  brick  elements  and  is  shown  in  Figure  5-1.  The  brick  elements  were  used  to 
provide  better  detail  in  modeling  the  barrels.  The  second  version,  called  the  simple 
barrel,  used  39  elements  with  8  cross  sections  and  is  shown  in  Figure  5-2.  The 
latter  version,  with  its  fewer  elements,  is  the  best  choice  for  providing  a  smoother 
animation  picture  and  keeping  the  number  of  degrees  of  freedom  within  a 
manageable  budget.  The  advanced  barrel  version  was  kept  for  comparison 
purposes  so  that  when  futiue  barrel  designs  are  tested  using  this  model  of  the  gun, 
the  most  complete  finite  element  modeling  technique  can  be  used. 

The  first  component  of  the  PHALANX  gun  to  be  received  at  the  Naval 
Postgraduate  School  was  a  single  barrel.  It  provided  an  opportunity  for  a 
validation  test  of  the  gun  model.  The  barrel  was  set  up  in  a  mechanical  engineering 
lab,  where  it  was  suspended  with  two  low-stiffiiess  cords  to  determine  its  normad 
modes  of  vibration  [Ref.  17].  This  method  is  called  the  "firee-firee  configuration."  If 
the  normal  modes  were  found  while  the  barrel  was  clamped  at  its  base  end,  for 
example,  the  modes  would  change  to  match  the  clamped  boxindary  conditions  and 
would  not  represent  properties  of  the  barrel  above.  In  the  firee-firee  configuration 
the  modes  of  vibration  of  the  real  gun  barrel  can  be  determined  for  use  in  the 
computer  models. 

Using  an  instrumented  hammer,  the  gun  barrel  was  tapped  near  its  base. 
The  resulting  vibrations  were  recorded  opposite  the  impact  point  with  a  single 
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barrel  made  of  brick  elements. 


barrel  made  of  beam  elements. 


transducer  or  accelerometer.  Appendix  A  shows  a  plot  of  the  recorded  response  for 
frequencies  up  to  3000  Hz.  The  excitation  force  due  to  the  gun  firing  at  a 
firequency  of  75  Hz  defines  the  important  fiequency  range  for  this  problem.  The 
normal  modes  of  vibration  of  the  two  versions  of  the  gun  barrel  were  calculated  in 
the  free-firee  configuration.  Here,  kinematic  degrees  of  fireedom  were  used  to 
simulate  the  restraints  provided  by  the  two  low  stiffiaess  cords  in  the  experimental 
set-up.  Figure  5-3  shows  the  first  mode  of  vibration  of  the  advanced  barrel  and 
Figure  5-4  shows  that  of  the  simple  barrel.  Table  5-1  lists  the  experimentally 
determined  modes  of  vibration  of  the  real  barrel  and  those  calculated  for  the 
computer  models.  Figure  5-5  graphically  illustrates  this  data. 


Figure  5-3.  The  first  normal  mode  of  vibration  of  the  advanced  barrel. 


TABLE  5-1.  COMPARISON  OF  ACTUAL  AND  MODELED  GUN  BARRELS 


Modes  of  Vibration 

Experimental  Barrel  Data 

Simple  Barrel  Data 

Advanced  Barrel  Data 

72.8 

71.6 

64.8 

217 

224 

201 

431 

455 

400 

694 

744 

652 

1040 

874 

809 

1428 

1111 

982 

1848 

1515 

1357 

2363 

1550 

1501 

2913 

2013 

1775 

Looking  at  Figure  5-5,  it  can  be  seen  that  there  is  good  correlation  with  both 
models  and  that  of  the  real  barrel  over  the  first  four  modes  of  vibration  or  up 
through  roughly  700  Hz.  According  to  general  finite  element  modeling  practice, 
this  agreement  of  the  first  four  modes  suggests  that  both  gun  models  are  validated. 
The  agreement  through  700  Hz  for  both  gun  models  is  considered  excellent.  It  is 
also  interesting  to  note  that  the  simple  barrel  data  matched  the  experimental 
barrel  data  more  closely  at  the  higher  frequencies  than  the  advanced  barrel.  This 
was  probably  due  to  the  fact  that  the  single  layer  of  brick  elements  used  in  the 
advanced  barrel  model  were  not  representing  beam  behavior  very  well.  As 


Figure  5-4.  The  first  normal  mode  of  vibration  of  the  simple  barrel. 

mentioned  in  the  previous  chapter,  more  bricks  are  required  to  reflect  bending  and 
shear  than  do  beam  elements.  Therefore,  in  the  case  of  the  advanced  barrel,  there 
should  at  least  be  two  layers  of  brick  elements  used  to  define  the  barrel’s  thickness. 

B.  MODES  OF  VIBRATION  OF  THE  COMPLETE  GUN  MODEL 

After  completion  of  the  PHALANX  gun  model,  its  static  properties  were 
validated  against  the  actual  gun.  The  next  step  was  to  calculate  its  normal  modes 
of  vibration.  Table  5-2  lists  the  first  21  calculated  modal  frequencies.  Upon 
observing  the  normal  modes,  it  became  apparent  that  the  first  mode,  which  had  a 
frequency  around  0.3  Hz,  was  the  normal  mode  of  vibration  of  rotation  about  the 
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Modes  of  Vibration 


Figure  5>5.  Comparison  of  the  modes  of  vibration  of  the  real  barrel  versus 
the  two  computer  models. 


lengthwise  axis.  This  is  a  computational  error  since  the  gun  barrel  assembly  is 
free  to  spin  in  the  model  and  does  not  have  any  torsional  restoring  force  associated 
with  this  motion.  Therefore,  it  has  a  nonlinear,  rigid  body  rotational  mode  of 
vibration  of  essentially  0.0  Hz.  The  fact  that  it  was  not  exactly  zero  was  probably 
due  to  the  bearing  springs  not  having  coincident  nodes  which  creates  an  unrealistic 
moment  arm  (see  Chapter  VII.  Suggestions  for  Future  Work).  As  mentioned  in 
the  previous  chapter,  the  normal  modes  that  are  being  calculated  are  those  without 
gyroscopics  or  an  unbalance  (centrifugal)  excitation.  Therefore,  the  normal  modes 
of  the  gun  were  resolved,  this  time  restricting  the  rotational  motion  of  the  rotor, 
which  removed  that  first  torsional  mode.  Hence,  the  first  20  flexible  modes  were 
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calculated.  The  periodic  force  due  to  the  firing  of  the  gun  is  at  a  frequency  of  75 
Hz  and  lies  between  modes  6  and  7.  A  complete  pictorial  record  of  all  the  modes 
used  in  the  dynamic  analysis  of  the  gun  appears  in  Appendix  B. 

It  is  also  important  to  note  that  in  the  n-^de  shape  pictures  to  follow,  the 
amplitude  of  the  deformation  is  greatly  exaggerated  in  order  to  the  general 

shape  of  the  normal  mode  and  the  relative  displacements.  Correct  amplitudes  will 
be  found  later  when  the  dynamic  response  of  the  system  is  evaluated  after  driving 
it  with  the  periodic  firing  force  and  values  of  damping  are  assigned.  Obtaining 
normal  modes  is  an  important  step  in  the  process  of  obtaining  displacement 
information. 


TABLE  5-2.  MODES  OF  VIBRATION  OF  THE  PHALANX  GUN 


Modes  of  Vibration 

Normal  Mode 

Frequency  (Hz) 

Normal  Mode 

Frequency  (Hz) 

1 

0 

12 

231 

2 

14.4 

13 

242 

3 

15.5 

14 

246 

4 

16.0 

15 

254 

5 

35.3 

16 

254 

6 

55.3 

17 

274 

7 

111 

18 

274 

8 

116 

19 

289 

9 

129 

20 

308 

10 

165 

21 

309 

11 

178 

1.  First  Flexible  Mode 

The  first  mode,  at  a  frequency  of  14.4  Hz,  is  considered  to  be  significant 
in  the  analysis  of  the  dispersion  data.  This  is  the  second  normal  mode  in  Table  5- 
2.  It  appears  to  be  a  link  between  the  periodic  excitation  force  and  its  rotation 
which  may  be  causing  a  resulting  whirling  motion  of  the  six-barrel  system.  Here, 
the  barrels  act  like  an  overhung  shaft  which  whips  around  its  center  line  with  the 
barrel  tips  tracing  a  new,  larger  circle.  This  ty]7e  of  motion  may  provide  an 
explanation  for  dispersion  patterns  that  appear  to  be  a  circle  of  the  six  individual 
barrel  patterns  as  mentioned  in  the  first  chapter.  This  motion  will  be  discussed  in 
greater  detail  in  the  last  chapter. 
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Figtire  5*6  shows  a  snapshot  of  this  first  flexible  mode.  It  primarily 
consists  of  a  rigid  body  motion  where  the  barrels  are  straight  and  pitch  from  side 
to  side.  The  rotor  through  the  barrels  move  with  respect  to  a  stationary  gun  body. 
These  components  are  pivoting  about  the  thrust  and  needle  bearings.  There  is 
little,  if  any,  fore  and  ail  motion.  There  is  no  recoil  motion  as  indicated  by  the 
recoil  adapter  springs  which  are  not  exercised.  In  this  mode  and  all  the  modes  to 
follow,  one  component  that  appears  to  be  "taking  a  lot  of  stress"  is  the  rear  gun 
body  wall.  Its  deformation  is  readily  apparent  when  a  transparent  picture  is  taken 
showing  its  motion  with  respect  to  the  rotor,  as  can  be  seen  in  Figure  5-7. 


Figure  5-6.  Hidden  line  image  of  the  first  flexible  mode. 
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Figure  5-7.  Hidden  line  image  showing  the  deformation  of  the  rear  gun 
body  plate. 

2.  Second  Flexible  Mode 

The  second  flexible  mode,  at  15.5  Hz,  is  a  combination  of  fore  and  aft 
motion  and  pitching  of  the  gun  barrel  assembly.  This  is  the  third  normal  mode  in 
Table  5-2.  Figure  5-8  shows  that  the  barrel  assembly  remains  rigid  as  it  pitches  in 
the  vertical  plane.  At  the  same  time,  the  whole  gun  recoils  and  exercises  the  recoil 
adapter  springs.  It  is  also  apparent  that  the  gun  body  is  also  being  exercised  in 
this  mode  for  it  shows  some  deformation. 

This  is  one  mode  where  an  easy  hand  calculation  can  be  done  in  order  to 
ensure  that  the  numerical  value  of  the  normal  mode  frequency  is  reasonable. 
Assuming  that  the  rotor,  barrel  assembly  and  gun  body  are  rigid,  the  motion  can 
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Figure  5-8.  Hidden  line  image  of  the  second  flexible  mode. 

be  described  as  a  symmetric  mode  of  the  recoil  adapters.  That  is.  the  whole  gun 
of  mass  M  =  120  kg  is  attached  to  ground  via  the  two  recoil  adapter  springs  as 
shown  in  Figure  5-9. 

The  equation  of  motion  for  such  a  system  is 

iWi' +2K)e{Xo -x)  =  0  .  (5.1) 

Here,  Xo  is  the  equilibrium  position  of  the  lumped  mass  before  it  vibrates  a 
distance  x  along  the  x  axis.  Kr  =  590,485  cN/cm,  is  the  stiilhess  value  for  the 
recoil  adapter  springs  used  in  the  numerical  model.  Prom  the  equation  of  motion, 
the  normal  mode  angular  fi*equency  of  vibration,  (o,  is 
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This  compares  favorably  with  the  model  value  of  f  =  15.5  Hz.  The  analytic 
calculation  assumed  a  mode  of  oscillation  purely  along  the  x  axis,  whereas  in  the 
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numerical  model,  there  is  a  pitching  motion  aa  well  as  the  recoil  motion  which 
would  contribute  to  the  determination  of  its  modal  frequency. 


3.  Third  Flexible  Mode 

The  third  flexible  mode,  or  fourth  normal  mode,  is  similar  to  the  previous 
mode.  A  snapshot  of  its  motion  is  shown  in  Figure  5-10.  Its  frequency  is  only 
slightly  greater  at  16.0  Hz.  The  one  apparent  difference  is  that  there  is  less  recoil 
motion  than  seen  in  the  previous  mode.  This  is  evidenced  by  the  lack  of 
deformation  in  the  recoil  adapter  springs.  Its  motion  is  almost  a  purely  pitching 
motion  of  the  barrel  assembly  in  the  vertical  direction. 


Figure  5-10.  Hidden  line  image  of  the  third  flexible  mode. 


4.  Fourth  Flexible  Mode 


The  first  thing  noticeably  different  in  the  fourth  flexible  mode,  or  fifth 
normal  mode,  is  the  massive  deformation  of  the  gun  body  as  displayed  in  Figure  5- 
11.  The  mode  of  oscillation  of  the  g\in  body  shows  a  bending  node  in  its  mid- 
section  with  a  normal  mode  frequency  of  35.3  Hz.  The  gun  barrel  assembly 
remains  rigid  with  respect  to  the  thrust  and  needle  bearings.  The  motion  of  the 
barrel  assembly  appears  to  be  caused  solely  by  the  bending  of  the  gun  body  as  it 
heaves  up  and  down.  There  is  no  recoil  motion  in  this  mode.  It  is  important  again 
to  note  the  deformation  of  the  rear  gun  body  plate  which  seems  to  suggest  that  its 
strength  would  play  a  crucial  role  in  restricting  barrel  tip  displacement. 


Figure  5-1 1.  Hidden  line  image  of  the  fourth  flexible  mode. 
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5.  Fifth  Flexible  Mode 


The  fifth  flexible  mode,  or  fourth  normal  mode  is  at  55.3  Hz.  Again  it 
appears  to  involve  a  mode  of  oscillation  of  the  gim  body  alone  as  shown  in  Figvire 
5-12.  However,  the  gun  body  oscillation  is  that  of  a  breathing  mode  where  the 
walls  move  in  and  out.  Like  the  previous  mode  the  gun  barrel  assembly  remains 
stationary  with  respect  to  the  bearings.  So  the  barrel  motion  is  directly  related  to 
the  motion  of  the  gun  body  walls.  As  the  walls  breath,  the  recoil  adapters  appear 
to  twist  in  the  vertical  direction.  There  is  no  motion  in  the  fore  and  aft  direction 
except  that  caused  by  gun  body  oscillation. 


Figure  5-12.  Hidden  line  image  of  the  fifth  flexible  mode. 
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6.  Sixth  Flexible  Mode 


The  final  mode  to  be  discussed  is  the  seventh  normal  mode  and  is  the 
first  mode  in  which  the  barrels  themselves  are  exercised.  This  mode,  at  1 1 1  Hz,  is 
shown  in  Figure  5*13.  The  barrel  tips  remain  parallel  to  ground  as  they  oscillate 
about  a  single  bending  node  located  just  forward  of  the  mid-barrel  clamp.  The 
boundary  conditions  of  the  muzzle  clamp  seems  to  be  vital  in  relation  to  barrel  tip 
motion.  In  the  real  gun,  three  of  the  barrels  are  unlocked  in  the  muzzle  clamp 
while  the  other  three  are  locked.  In  this  model,  all  the  barrels  are  rigidly  attached 
or  locked  to  the  muzzle  clamps.  A  possible  future  modification  would  be  to  unlock 
three  of  the  barrels  in  the  muzzle  clamp  and  simulate  the  effect  of  a  gap  around 
them  using  spring  elements.  The  difficulty  would  come  in  determining  the  stiffness 
values  of  the  springs  since  this  is  a  nonlinear  problem.  A  new  version  of  I-DEAS, 
called  the  Master  Series,  is  able  to  simulate  the  effect  of  a  gap  or  tolerance.  After 
the  real  g\in  is  set  up  on-sight,  then  tests  can  be  performed  to  determine  the 
proper  stiffness  values.  This  vibrational  mode  of  the  barrels  with  damps  seems  to 
exhibit  the  behavior  of  a  deep  flange  beam  which  is  different  than  the  six  barrels 
represented  as  a  single  beam.  A  deep  flange  beam  tends  to  shear  like  the  barrels 
in  the  model  rather  than  bend  like  a  beam. 

Looking  at  the  rest  of  the  gun,  there  is  massive  deformation  of  the  gun 
body  with  a  single  bending  mode  in  the  middle.  The  rotor  essentially  remains 
parallel  to  ground  as  the  gun  body  and  barrels  flex.  There  is  no  recoil  motion 
assodated  with  this  normal  mode. 

C.  MODAL  TRUNCATION 

With  22,470  degrees  of  freedom,  there  are  an  equal  number  of  normal  modes 
of  vibration.  Solving  for  all  of  these  normal  modes  is  unnecessary,  but  would 
defeat  the  purpose  of  reducing  the  complex  model  into  relatively  simple  normal 
modes.  The  System  Dynamics  Analysis  module,  to  be  discussed  fiirther  in  the  next 
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Figure  5*13.  Hidden  line  image  of  the  sixth  flexible  mode. 

chapter,  calculates  the  dynamic  response  of  a  system  by  first  solving  for  the  modes 
of  vibration.  Instead  of  using  the  22,470  physical  degrees  of  freedom,  it  can 
describe  the  response  using  the  first  20  flexible  modes  of  Table  5-2. 

According  to  generally  accepted  finite  element  modeling  practice,  the  number 
of  normal  modes  that  should  be  carried  through  for  a  dynamic  response  analysis 
are  those  that  fall  within  twice  the  excitation  firequency.  Since  the  gun’s  principal 
excitation  fiiequency  is  75  Hz,  the  modes  through  150  Hz  should  be  calculated.  By 
solving  for  the  first  20  modes,  up  to  a  frequency  of  309  Hz  then  this  criteria  i 
more  than  satisfied. 


87 


VI.  DRIVING  AND  DAMPING  OF  THE  PHALANX  GUN 


A.  DYNAMIC  ANALYSIS  OF  THE  GUN  MODEL 

Having  completed  a  modal  analysis  of  the  PHALANX  gun  model,  the  next 
step,  as  outlined  in  Chapter  III,  is  to  perform  a  forced  response  analysis  of  the 
gun.  The  excitation  forces  caused  by  the  bullets  recoil  and  estimates  of  damping  in 
the  gun  are  applied.  The  resulting  response,  such  as  barrel  tip  displacement,  can 
then  be  determined  and  compared  to  experimental  data.  The  most  prominent 
normal  mode  can  be  found  to  describe  the  actual  motion  of  the  gun  when  firing. 
With  the  completion  of  these  steps,  a  working  model  of  the  PHALANX  gun  has 
been  created.  The  factors  which  contribute  to  dispersion  can  be  studied 
systematically,  and  potential  solutions  to  problems  can  be  tested  to  provide  an 
additional  source  of  information  to  the  PHALANX  program  office.  Various 
components  can  be  "swapped"  out  to  evaluate  modified  components,  and  the 
resulting  effects  on  gun  performance  predicted  prior  to  any  hardware  tests.  This  is 
the  first  finite  element  model  of  its  kind  ever  constructed  of  the  PHALANX  gun, 
and  the  possibilities  for  its  use  are  only  limited  by  imagination. 

The  I-DEAS  computer  code  evaluates  a  structxire’s  response  to  applied  forces 
in  the  computer  module  called  System  Dynamics  Analysis  (SYSTAN).  The  method 
used  by  SYSTAN  is  to  first  develop  a  modal  component  model  of  the  structure. 
That  is,  a  certain  number  of  normal  modes  of  a  structure  are  calculated  and  stored 
as  a  new  modal  model  of  the  structure  which  is  then  used  to  calculate  the  dynamic 
response.  Other  methods  can  be  used,  such  as  directly  integrating  the  equation  of 
motion 


F  = 


M 


A'  + 


X  + 


K 


(6.1) 


Here,  column  matrices  are  represented  with  underlines  and  square  matrices  with 
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brackets.  F  represents  the  applied  forces  and  their  direction  on  each  degree  of 
freedom  in  the  system.  [M]  is  the  modal  mass  matrix  that  was  obtained  through 
the  process  of  transforming  to  normal  coordinates  described  in  Chapter  III.  Each 
of  the  ith  matrix  elements  is  a  measure  of  how  much  mass  is  participating  in  the 
ith  normal  mode.  Similarly,  [K]  is  the  modal  stiffness  matrix  and  is  a  measure  of 
the  stiffness  occurring  in  each  of  the  normal  modes.  [C]  is  the  associated  damping 
for  and  between  each  degree  of  freedom.  The  displacements,  velocities  and 
accelerations  for  each  degree  of  freedom  are  represented  by  X  ,X  ,  and  X 
respectively.  However,  when  there  is  a  structure  with  many  degrees  of  freedom, 
like  the  PHALANX  gun  model  with  22,470,  it  becomes  apparent  that  the 
calculations  required  are  memory  and  time  consuming.  Reducing  the  number  of 
physical  degrees  of  freedom,  X,  &x>m  22,470  to  20  modal  degrees  of  freedom,  x, 
greatly  reduces  the  complexity  of  the  model.  In  order  to  accomplish  this  reduction, 
a  transformation  to  normal  coordinates  diagonalizes  the  mass,  stiffness  and 
damping  matrices  of  the  equation  of  motion  as  discussed  in  Chapter  III. 
Reductions  which  do  not  diagonalize  are  also  possible,  but  will  not  be  discussed 
here.  The  equation  of  motion  can  then  be  written  as 


o 

_ 1 

Cii  0  . 

Kn  0  . 

0  ^22  0 

0  ^22  0 

0  K22  0 

0  .  0 
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0 

.  0  My 

•  0  Cy 
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This  equation  is  similar  to  the  diagonalized  equation  of  motion  (3.60),  developed  in 
Chapter  III,  except  that  now  a  damping  matrix  [C  ]  and  an  excitation  force  F  have 
been  added.  Damping  is  necessary  when  evaluating  a  dynamic  response  in  order 
to  simulate  the  effects  of  energy  dissipation  in  each  normal  mode.  Otherwise,  the 
motion  of  a  structure  will  continue  to  grow  unbounded  as  the  excitation  force  is 
applied.  In  (6.2)  and  in  SYSTAN,  the  damping  is  assumed  to  be  viscous  damping 
proportional  to  the  velocity  of  the  structure  simulating  the  effects  of  friction.  With 
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the  matrices  diagonalized,  and  hence  uncoupled,  the  equations  of  motion  are  easily 
solved.  [Ref.  11] 


The  conversion  can  then  be  made  back  to  the  physical  deg^rees  of  freedom 
through  the  relation 


^  = 


(6.3) 


where  [<1>]  represents  the  spatial  and  temporal  parts  of  the  mode  shape  as 
discussed  in  Chapter  III.  Each  column  of  this  matrix  is  an  eigenvector  0‘  which 
corresponds  to  each  modal  degree  of  freedom.  [Ref.  11] 


B.  DAMPING  IN  THE  PHALANX  GUN 

Damping  in  the  PHALANX  gun  model  is  defined  as  Rayleigh  or  proportional 
damping.  Since  the  diagonalized  modal  mass  and  stiilhess  matrices  were  obtained 
without  the  velocity  damping  term  iu  the  equation  of  motion,  then  the  damping 
matrix  [C]  is  not  diagonahzed.  Instead,  the  damping  coefficients  are  defined  by 
the  relation 

C=oJir  +  pAf  (6.4) 

Diagonalizing  this  matrix  yields 

0^C0  =  +  p0^Af0  (6.5) 

=  a[K]  +  p[Ml  , 

where  [K]  and  [M]  are  now  diagonalized  matrices.  Substituting  these  relations 
back  into  the  equation  of  motion  gives 

F  =  [M]i  +  [K]x  +  (a[K]  +  p[M])x  .  (6.6) 

The  quantity  in  parenthesis  is  called  modal  damping  and  it  was  this  quantity  that 
was  assigned  a  nximerical  value  in  the  gun  model.  The  correct  modal  damping 
values,  a  and  P,  can  only  be  determined  by  measurements  on  the  actual  system 
being  modeled.  Therefore,  damping  is  associated  with  each  normal  mode  and  not 
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with  each  element  in  the  gun  model.  To  more  accurately  represent  damping  of  the 
real  gun,  damping  would  have  to  be  defined  in  terms  of  two  physical  points,  such 
as  a  dashpot  in  the  recoil  adapter,  and  not  just  with  each  mode.  Each  mode  itself 
is  not  a  pure  motion  as  well.  For  example,  the  recoil  mode  has  a  certain  amount  of 
pitching  associated  with  it,  so  an  additional  value  of  damping  should  be  assigned  to 
this  additional  motion.  In  the  absence  of  measurements,  the  particular 
components  of  the  gun  involved  in  each  normal  mode  can  be  iised  to  estimate 
damping.  Damping  values  were  assigned  based  on  knowledge  of  how  the 
components  are  assembled,  for  example,  welded  or  bolted  components  have  certain 
damping,  while  the  barrels  alone  have  only  small  damping  as  observed  in  the 
modal  vibration  tests  and  the  recoil  adapter  springs  are  intentionally  damped  by 
the  manufacturer,  as  described  by  a  hysteresis  plot. 

In  the  SYSTAN  module,  damping  is  defined  as  a  percent  of  critical  damping. 
To  clarify  further,  the  effects  of  damping  are  associated  with  three  physically 
distinct  cases.  These  cases  are  derived  firom  the  equation  of  motion  (6.2)  where  an 
exponential  function  of  the  form 

X=e’‘  ,  (6.7) 


was  used  as  a  solution  with  the  driving  force  set  to  zero.  The  equation  of  motion 
for  the  first  mode  is  then 


+  Kue’*  =  0  .  (6.8) 

dt^  dt 

Performing  the  necessary  derivatives  and  canceling  the  factors  of  yields  the 
quadratic  equation 

+  Ciig  +  Kii  =  0  (6.9) 


Solving  for  q  using  the  quadratic  formula  gives 

X 

_  -C|i  ±  -  4M^iKn)^ 

■  2Mh 


(6.10) 
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The  three  difTerent  damping  cases  are  then  determined  from  the  value  of  the 
discriminant  Cn^  -  4MiiKii.  Ouerdamping  occurs  for  the  case  where 
>  0.  In  this  case  there  is  an  exponential  type  of  decay.  Given  an 
initial  displacement,  the  structure  does  not  oscillate  about  its  equilibrium  position 
as  shown  in  Figiire  6- Kb).  The  critical  damped  case  occurs  when 

Cii^  -  =  0.  Like  overdamping,  critical  damping  does  not  oscillate,  but 

dampens  out  asymptotically  as  shown  in  Figxire  6-l(c).  Finally,  the  underdamped 
case  occurs  when  Cn^  -  <  0  where  the  structure  oscillates  about  its 

equilibrium  point  sinusoidally,  but  eventually  dampens  out  as  shown  in  Figirre  6- 
1(a).  [Ref.  18] 


THE  DAMPED  HARMONIC  OSOULATOR 


t 


Remm  of  hamioaic  oKillitor  to  eqoUibnnin.  (a)  Underdampcd.  (b)  Ovudampod. 
|c)  CritiaU;  damped. 


Figure  6-1.  Graphs  of  displacement  for  the  three  case  of  damping. 
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When  determining  a  damping  value  based  on  knowledge  of  how  the 
components  are  connected,  experts  in  the  industry  have  accepted  general  empirical 
values  until  more  information  is  obtained.  Modes  involving  welded  parts  are 
assumed  to  have  damping  values  of  around  1  percent  of  critical  damping  whereas 
bolted  structures  have  values  between  3  and  5  percent.  For  clamped  components, 
2  percent  is  an  accepted  value.  Since  most  of  the  components  are  clamped  in  the 
gun,  2  percent  was  used  for  most  of  the  modes.  [Ref.  11] 

A  second  case  estimates  the  modal  damping  value  for  a  single  gun  barrel. 
Only  this  component  has  been  observed  while  measuring  the  modal  vibration. 
While  obtaining  the  normal  modes  of  vibration  of  the  single  g\jn  barrel  by  tapping 
it  with  the  instrumented  hammer,  it  was  observed  to  have  a  very  low  damping 
value.  It  was  observed  to  have  a  damping  value  of  less  than  1  percent  in  the  free- 
free  configuration.  The  damping  value  used  in  the  computer  model  was  raised  to  2 
percent  to  reflect  the  fact  that  each  gun  barrel  is  clamped  at  its  base,  middle  and 
end  points. 

Another  component  that  has  experimental  damping  information  available  are 
the  two  recoil  adapters.  A  hysteresis  plot  provided  by  Martin  Marietta  for  the  two 
recoil  adapters  is  shown  in  Figure  6-2  [Ref.  15].  It  shows  the  loading  force  applied 
to  the  recoil  springs  for  one  cycle  on  the  vertical  axis  and  the  resulting  deflection  on 
the  horizontal  axis.  The  upper  branch  reflects  the  initial  loading  and  the  lower 
branch  the  return  to  battery.  The  area  contained  within  the  two  branches  is  the 
energy  dissipated  per  cycle,  £/)=  4.34  x  10^ N-cm  (384  lb -in).  The  damping 
coefficient ,  Cfty,t ,  is  then  determined  from  the  equation 

Ed  =  ,  (6.11) 

where,  x  is  the  amplitude  of  oscillation  about  its  driven  equilibrium  point 
estimated  by  the  mid- point  in  the  hysteresis  plot,  or  x  =  0.4  cm  (0.15  in).  After  the 
first  few  rounds  are  fired,  the  gun  sets  back  about  1.3  cm  (0.5  in),  and  then 
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oscillates  about  this  new  equilibrium  position  xintU  the  gun  stops  firing.  It  never 
returns  to  battery  between  shots  Uke  other  naval  guns.  The  angular  frequency  of 
oscillation  of  the  recoil  mode  O)  =  30it/  s .  Rearranging  and  solving  for  yields 


C/tyat 


=  91.6 

jj Lr  ro)  cm!  s 


(6.12) 


The  damping  coefficient,  y=  C}^y^l  Ccru>  used  by  the  computer  code  is  expressed  as 
a  percentage  of  critical  damping,  Cent  where 

=  AmK  ,  (6.13) 


m  is  the  mass  of  the  gun,  and  K  is  the  stiffness  of  the  recoil  adapters. 
Substituting  values  into  this  equation  gives 


=  2.38  X  10^- 


N 


cm  I  s 


(6.14) 


so  that  y  can  now  be  calculated  as 


Y=^  =  0.39 

^erit 


(6.15) 


The  value  for  y  in  (6.15)  is  consistent  with  the  actual  damping  value  used  in 
the  design  of  the  recoil  adapters,  y  =  0.33,  as  provided  by  Martin  Marietta  [Ref. 
19].  A  damping  coefficient  which  approaches  critical  damping,  y  1,  gives  a  more 
optimum  return  to  equilibrium  after  the  excitation  force  has  stopped.  In  the 
computer  model  a  damping  value  of  y=  0.33  was  used  for  this  mode.  The  recoil 
mode,  referred  to  in  Chapter  V  as  the  third  flexible  mode  of  vibration  at  15.5  Hz, 
shows  the  gun  sliding  fore  and  aft  exercising  the  recoil  adapter  springs  with  a 
slight  pitching  of  the  barrels  in  the  vertical  plane.  This  is  the  only  mode  using  a 
damping  value  other  than  2  percent.  The  purpose  of  the  recoil  adapters  is  to 
dissipate  some  of  the  force  caused  by  the  firing  rounds  and  to  couple  the  remaining 
force  to  the  ship.  Allowing  a  gun  to  recoil  even  a  few  calibers  greatly  reduces  the 
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DEFLECTION  IN  INCHES 

Figure  6-2.  Hysteresis  plot  of  the  recoil  adapter  springs. 

force  of  recoil  transmitted  to  the  mount  so  that  a  much  Ughter  gun  mount  can  be 
used. 

Table  6-1  lists  the  damping  values  applied  to  each  normal  mode  of  vibration 
used  in  the  PHALANX  gun  model.  Most  of  the  damping  values  used  in  the  model 
are  crude  estimates  based  on  information  derived  without  the  actual  gun  on  sight. 
Further  experimental  tests  should  be  performed  allowing  for  further  refinement  of 
the  damping  values. 
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TABLE  6-1.  MODAL  DAMPING  VALUES 


Damping  of  the  Modes  of  Vibration 

Normal  Mode 

Damping  (y) 

Normal  Mode 

Damping  (y) 

1 

0 

12 

0.02 

2  (yaw) 

0.02 

13 

0.02 

j 

3  (recoil,  pitch) 

0.33 

14 

0.02 

4  (pitch) 

0.02 

15 

0.02 

5  (gun  body  mode) 

0.02 

16 

0.02 

6  (gun  body  mode) 

0.02 

17 

0.02 

7  (barrel  bending) 

0.02 

18 

0.02 

8 

0.02 

19 

0.02 

9 

0.02 

20 

0.02 

10 

0.02 

21 

0.02 

11 

0.02 

C.  DRIVING  THE  PHALANX  GUN 

The  principal  excitation  force  of  the  PHALANX  gun  is  the  repeated  firing  of 
the  rounds.  Only  one  barrel  shoots  at  a  time,  so  that  a  firing  rate  of  4500  shots 
per  minute  gives  the  frequency  of  the  excitation  force  75  Hz.  The  experimentally 
determined  force  on  the  breech  for  a  single  shot  lasts  roughly  2.3  milliseconds  as 
displayed  in  Table  6-2  and  plotted  in  Figure  6-3  [Ref.  15]. 
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TABLE  6-2.  FORCE  ON  GUN  BREECH 


Force  of  Single  Shot  on  Breech 

Time  (milhseconds) 

Force  (lb) 

Time  (milhsecond) 

Force  (lb) 

0.05 

8.74  X  10^ 

0.80 

2.85  X  10“ 

0.10 

1.26  X  10^ 

0.85 

3.02  X  10“ 

0.15 

1.72  X  10^ 

0.90 

3.12  X  10“ 

0.20 

2.33  X  10^ 

0.95 

3.16  X  10“ 

0.25 

3.12  X  10^ 

1.00 

3.13  X  10“ 

0.30 

4.14  X  10^ 

1.05 

3.02  X  10“ 

0.35 

5.45  X  10^ 

1.15 

2.59  X  10“ 

0.40 

7.11  X  10^ 

1.25 

2.02  X  10“ 

0.45 

9.17  X  10^ 

1.45 

1.16  X  10“ 

0.50 

1.17  X  10“ 

1.65 

7.68  X  10^ 

0.55 

1.45  X  10“ 

1.85 

5.57  X  10^ 

0.60 

1.76  X  10“ 

2.05 

4.29  X  10^ 

0.65 

2.08  X  10“ 

2.25 

3.46  X  10^ 

0.70 

2.37  X  10“ 

2.31 

3.26  X  10^ 

0.75 

2.64  X  10“ 

A  t3rpical  missile  engagement  burst  lasts  for  only  a  few  seconds  up  to  a 
maximum  of  5  seconds  before  the  barrels  deteriorate  due  to  the  high  temperature. 
Therefore,  a  driving  force  of  one  hundred  shots  was  constructed.  Due  to  the 
sharpness  of  the  firing  "spike"  compared  to  the  length  of  time  between  shots,  13.3 
milliseconds  at  75  Hz,  it  was  necessary  to  determine  more  data  points  for  each 
shot  in  order  to  increase  the  resolution  of  the  response.  A  function  was  fit  to  the 
available  data  in  Figure  6-3,  as  shown  in  Figure  6-4.  The  function’s  representation 
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Figure  6-3.  Plot  of  the  force  on  the  breech  from  a  single  shot  [Ref.  15). 
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of  the  first  3  shots  is  shown  in  Figure  6-5  with  dots  representing  the  experimental 
data  points  described  above.  The  function  was  used  in  SYSTAN  to  create  a  forcing 
function  using  26,400  data  points  over  a  one  hundred  shot  burst.  The  26,400 
evenly  spaced  points  are  spaced  at  intervals  of  0.05  milliseconds  to  provide  a  good 
representation  of  each  2.3  millisecond  shot.  The  number  of  data  points  used  is 
then  determined  as 


Number  of  points 


- ^-3  seconds  _  26,400 

0.05  milliseconds 


where,  1.3  seconds  is  the  length  of  time  for  a  100  shot  burst. 


(6.16) 
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Figure  6-4.  Plot  of  the  function  used  to  approximate  the  single  shot  data. 

A  single  shot  hy  itself  is  broad  band  and  is  comprised  of  many  frequencies.  It 
is  Uke  an  impulse  or  delta  fimction.  The  firing  of  a  single  shot  by  the  gun  model, 
would  be  expected  to  excite  many  normal  modes  of  vibration  over  a  wide  frequency 
range.  However,  the  gun  is  never  used  in  the  single  shot  mode,  so  this  case  is 
irrelevant.  A  long  series  of  shots  is  no  longer  broad  band,  but  has  a  fimdamental 
frequency  of  75  Hz.  This  type  of  excitation  force  would  be  expected  to  excite  fewer 
modes  in  a  small  range  around  75  Hz. 

A  burst  of  shots  is  applied  at  the  base  of  the  firing  barrel  located 
approximately  between  10  and  12  degrees  from  top  dead  center.  The  firing 
sequence  was  first  designed  for  the  older  firing  rate  of  3000  shots  per  minute.  In 
this  case,  the  round  was  ignited  at  this  off-center  position  so  that  when  the  barrel 
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Figure  6-5.  The  function’s  approximation  for  the  hrst  3  shots. 


reached  the  top,  the  round  would  just  be  leaving  the  barrel.  However,  the  firing 
rate  has  been  changed  to  4500  shots  per  minute  without  an  adjustment  being 
made  to  the  igniting  point.  Because  of  the  higher  rotation  rate  of  the  barrel 
assembly,  the  bullet  leaves  the  barrel  approximately  14  degrees  past  top  dead 
center.  The  time  varying  magnitude  of  the  firing  force  is  applied  to  the  gun  over  a 
20  degree  arc.  To  simulate  this  off-center  firing  force,  the  barrels  were  moved  10 
degrees  from  the  top.  The  gun  model  barrels  can  be  rotated  to  any  point  along  this 
20  degree  arc,  and  the  fiiring  force  reapplied.  The  barrel  tip  motion  can  then  be 
recalculated.  This  is  an  approach  being  taken  by  Martin  Marietta.  It  has  not  been 
possible  to  continuously  rotate  the  barrels  in  the  numerical  model  while  applying 
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the  firing  force.  The  current  arrangement  gives  a  good  approximation.  It  appears 
that  there  is  a  way  to  "trick”  the  computer  code  into  simulating  the  effects  of 
rotation,  but  this  is  ctirrently  beyond  the  scope  of  the  current  work.  Exploring  this 
improvement  is  recommended  as  a  further  refinement  of  the  model. 

A  possible  comphcation  in  applying  the  firing  force  to  the  gun  is  the  fact  that 
it  may  be  a  follower  force  problem.  That  is,  as  the  barrel  assembly  becomes 
displaced  from  its  center  line,  for  whatever  reason,  the  firing  force  ought  to  be 
applied  to  the  gun  in  this  displaced  position.  The  combination  of  this  displacement 
and  normal  firing  would  add  other,  currently  unpredictable,  effects  on  the  barrel 
tip  motion.  The  following  force  problem  is  probably  best  suited  to  be  described 
analytically.  Preliminary  steps  in  this  direction  have  been  taken  by  Professor 
Gordis. 

Another  force  acting  on  the  gun  arises  firom  the  rotating  barrels.  The  rotation 
rate  of  the  six-barrel  system  is  12.5  Hz  =  75  Hz/  6  barrels.  This  motion  can  play 
an  important  role  in  gun  dynamics  since  there  is  a  nearly  resonant  normal  mode  of 
vibration  at  14.4  Hz.  This  normal  mode  describes  a  pitching  motion  of  the  barrels 
and  may  contribute  to  the  dispersion  pattern  if  excited. 

When  the  rotating  barrels  are  displaced  firom  their  center  line,  they  will  cause 
excitation.  A  perturbing  force  that  could  cause  a  displacement  would  be  a  weight 
imbalance  between  the  barrels.  The  mass  of  a  barrel  may  vary  by  as  much  as 
0.227  kg  (0.5  lb)  away  from  the  average  weight  of  8  kg  (18  lb).  This  imbalance 
would  cause  the  barrel  cluster  to  bend  slightly  as  they  rotate.  This  displacement 
could  also  start  when  the  barrel  cluster  kicks  up  as  the  gun  first  fired.  The  kicking 
up  is  caused  by  the  placement  of  the  recoil  adapters  and  the  rear  ball  joint  in 
relation  to  the  point  of  application  of  the  firing  force.  Displacement  is  created  due 
to  the  imbalance  of  the  firing  point  and  the  three  supports.  Any  of  these  causes 
leads  to  excitation  of  a  physical  motion  called  gyroscopic  "whirl."  The  near 
resonance  of  the  barrel  rotation  rate  and  the  first  flexible  mode  of  vibration 
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appears  sufficient  to  excite  this  pitching  mode  and  greatly  edTect  dispersion. 
Furthermore,  a  non-preloaded  thrust  bearing,  as  is  used  in  the  PHALANX  gun,  is 
considered  to  be  a  root  cause  of  gyroscopic  whirl  [Ref.  20].  Gyroscopic  whirl  will  be 
further  discussed  in  a  latter  section.  There  are  a  variety  of  "causes  '  for  dispersion 
that  may  feed  upon  others.  If  one  could  determine  the  primary  causes  then  it  may 
be  possible  to  prevent  or  limit  others. 

Another  excitation  force  on  the  gun  is  the  torque  exerted  on  the  barrel  as 
each  round  spins  up  as  it  travels  along  the  barrel.  The  barrel  lands  have  a  right- 
hand  exponential  gain  twist  which  rotates  the  penetrator  to  a  rotation  rate  of 
120,000  revolutions  per  minute  (2,000  revolutions/  second)  as  it  leaves  the  barrel 
[Ref.  21].  The  torque  to  achieve  this  rotation  rate  can  be  calculated  by  first 
determining  its  rate  of  angular  acceleration. 


t. 


(6.17) 


where,  (O/  is  the  final  rotation  rate,  ©o  is  the  initial  rotation  rate,  and  is  the 
time  that  the  round  travels  the  length  of  the  barrel.  Substituting  the  values  in 
yields 

i,  =  5.46  X  10' . 

second^ 

The  torque  applied  by  the  barrel  on  the  round  to  achieve  this  rotation  rate  can  be 
calculated  using  the  following  relation. 


(6.18) 


where  I  is  the  moment  of  inertia  of  a  cylinder  which  approximates  that  of  the 
roxind,  and  a  is  the  inner  radius  of  the  barrel.  I  is  determined  according  to  the 
relation 


rrira^ 

2 


(6.19) 
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where  is  the  mass  of  the  round.  /  is  calculated  to  be 

,  ^  (0.071  tgxio  mmf  ^  3  x  10-»  , 

A 

Substituting  this  value  into  equation  (6.18)  gives 

Fft  =  435.6  lb  , 

showing  a  significant  force  on  the  barrels.  This  type  of  force  may  be  difficult  to 
simulate  with  this  model  due  to  its  variation  in  location  with  time,  and  is  similar  to 
the  rotating  barrel  problem. 

Other  possible  forces  include  an  imbalanced  load  on  the  gun  with  heavy 
preshot  rounds  on  one  side  and  the  li|^ter  expended  casings  on  the  other.  Also, 
the  act  of  loading  a  round  is  not  insignificant  and  would  add  an  impulsing  torque 
to  the  gun  assembly.  Furthermore,  the  actual  rate  of  firing  varies  since  the 
hydraulic  motor  which  rotates  the  harrel  assembly,  and  determines  the  firing  rate, 
varies  its  speed  by  as  much  as  ±  50  revolutions  per  minute.  These  variations  could 
occur  as  sudden  rotation  speed  shifts  resulting  in  a  torque  or  "jerk"  on  the  barrel 
assembly. 

The  modeling  of  the  PHALANX  gun  in  this  study  ends  at  the  three 
supporting  fixtures.  Other  factors  affecting  the  gun’s  motion  from  beyond  these 
points  are  not  being  simulated,  but  may  also  play  a  role  in  the  actual  dispersion. 
The  gun  mount  will  rarely  be  in  a  stationary  position  when  attempting  to  shoot 
down  targets  moving  at  speeds  on  the  order  of  a  few  Mach.  The  gun  can  slew  at  a 
maximum  rate  of  92  degrees  per  second  in  elevation  and  126  degrees  per  second  in 
train  or  rotation.  This  motion,  may  serve  as  a  perturbing  force  to  start  the 
gyroscopic  whirl  mentioned  earlier.  The  whole  gun  mount  is  mounted  on  a  set  of 
shock  absorbers  whose  purpose  is  to  limit  the  shock  wave  of  a  close-order 
explosion.  The  vibrating  gun  interacting  with  these  shock  mounts  may  merit 
future  consideration. 
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Ib  an  actual  gxin,  all  of  these  forces  may  play  a  role  in  the  observed  dispersion 
pattern.  Some  of  these  forces  may  be  added  to  the  current  gun  model  analysis,  but 
this  study  only  addresses  the  role  played  by  the  most  significant  force  of  the  gun 
firing. 

D.  RESPONSE  OF  THE  GUN  MODEL 

While  driving  the  gun  model  with  a  one  hundred  shot  forcing  function,  the 
resulting  response  of  the  top  barrel  tip  was  evaluated  to  its  translation  along  each 
of  the  three  coordinate  axes.  Shown  in  Figure  6-6  is  the  amplitude  of  the  barrel 
tip’s  motion  along  the  x-axis  in  the  fore-aft  direction  versus  time.  The  plot  shows 
that  the  gun  immediately  sets  back  by  roughly  1.27  cm  (0.5  in)  after  the  first  shot 
is  fired.  It  then  begins  to  oscillate  about  its  new  eqiiilibrium  position  by  roughly 
0.80  cm  (0.32  in)  at  a  firequency  between  15  to  20  Hz  as  one  or  more  of  the  lower 
order  modes  is  excited.  It  appears  to  reach  a  steady-state  oscillation  near  the  end 
of  the  burst  of  shots.  When  the  firing  is  complete  after  1.3  seconds  time  on  the 
horizontal  axis,  the  gun  moves  forward  to  its  initial  position  and  vibrates  like  an 
underdamp)ed  harmonic  oscillator.  During  excitation,  the  displacement  versus  time 
shows  an  irregular  pattern  of  jagged  peaks  indicating  excitation  of  several  normal 
modes.  Understanding  which  normal  modes  are  involved  in  the  response  will 
provide  the  greatest  benefit  in  predicting  the  dispersion  pattern.  The 
determination  of  these  modes  will  be  provided  by  a  frequency  analysis  which  will 
be  the  subject  of  the  next  section. 

The  important  barrel  tip  displacement  in  the  vertical  direction  is  shown  in 
Figure  6-7.  The  barrel  tip  can  be  seen  to  pitch  upward  significantly  after  the  first 
shot,  and  then  approach  steady-state  oscillation  near  the  end  of  the  firing  burst. 
The  peak-to-trough  barrel  tip  displacement  measured  near  the  steady-state  motion 
is  about  2.5  mm  (0.1  in)  and  gives  the  observed  1.64  milliradian  dispersion  for  a 
1.52  m  long  barrel.  The  first  20  to  30  shots  of  the  firing  burst  have  a  much  greater 
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Figure  6-6.  Response  plot  of  the  barrel  tip  in  the  fore-aft  direction 
(x-plane). 

dispersion  pattern  than  the  last  20  to  30  shots.  No  data  is  currently  available  for 
comparison  to  this  modeled  phenomenon.  Since  the  recoil  normal  mode  of 
vibration  at  15.5  Hz  has  an  associated  pitching  motion,  it  appears  that  this  motion 
is  important  for  predicting  the  dispersion  pattern.  Unlike  the  response  in  the 
fore-aft  direction,  this  vertical  response  appears  to  involve  only  a  single  mode  with 
a  h  quency  of  oscillation  between  15  to  20  Hz.  The  barrel  tip  displacement 
corresponds  to  an  angle  of  1.64  milhradian  movement  which  represents  a 
significant  start  in  quantifying  the  amount  of  dispersion  associated  with  barrel  tip 
motion. 
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Figure  6-7.  Response  plot  of  the  barrel  tip  in  the  vertical  direction 
(y-plane). 

Figure  6-8  shows  barrel  tip  translation  along  the  z-axis,  or  side-to-side 
motion.  The  amplitude  of  motion  is  noticeably  less  than  in  the  vertical  plane  and 
demonstrates  the  relative  strengths  of  the  moment  arms  about  the  center  of 
gravity.  In  the  vertical  plane,  there  is  a  larger  moment  arm  from  the  firing  point 
to  the  recoil  adapter  pivot  points  than  there  is  in  the  z-plane  so  that  greater 
displacement  in  the  y-p^.-ne  is  expected.  The  off-center  firing  force  causes  a  larger 
displacement  in  one  direction  than  in  the  other.  This  shows  that  actual  dispersion 
should  be  broken  down  into  components  in  the  two  transverse  dimensions  x  and  z. 
Based  on  this  analysis,  the  Naval  Surface  Warfare  Center,  Dahlgren,  and  Martin 
Marietta  have  now  begun  to  describe  their  observed  dispersion  patterns  as 
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elliptically  shaped  [Ref.  22].  The  maximum  barrel  tip  displacement  in  the  z 
direction  is  about  5.10  mm  (0.20  in),  and  settles  into  a  steady-state  displacement 
of  about  2  mm  (0.08  in).  The  response  in  the  z-plane  appears  to  be  multimode 
which  is  similar  to  that  of  the  x-plane  response.  The  response  and  ring  down  show 
beats  which  suggests  that  the  multimode  response  involves  two  normal  modes 
whose  frequencies  are  close  together.  Recall  that  the  response  of  a  structure  is  just 
the  summation  of  the  excited  normal  modes  of  vibration  as  given  by 


X  = 


(6.20) 


The  beating  appearing  in  this  response  plot  indicates  that  the  normal  modes 
excited  have  similar  damping.  The  multiple  normal  modes  excited  in  the  x-plane 
apparently  had  different  damping.  More  will  be  discovered  in  the  next  section  on 
this  topic  using  the  firequency  response  plots. 


E.  FREQUENCY  RESPONSE  FUNCTIONS 

The  previously  discussed  response  plots  are  important  for  ensuring  that  the 
hnite  element  model  appears  to  be  a  close  approximation  of  the  real  structure.  It 
is  a  check  on  the  assvimptions  made  during  the  modeling  process.  However,  it  is 
more  important  to  understand  the  normal  modes  involved  in  the  response.  In  this 
way,  it  is  possible  to  pinpoint  causes  of  large  amplitude  displacements  of  the 
structure.  A  structure’s  response  to  an  excitation  depends  on  the  magnitude  of  the 
excitation  force,  the  frequency  of  the  excitation  force,  the  point  of  application  and 
damping.  Large  displacements  can  occur  if  the  firequency  of  the  excitation  is  near 
that  of  a  normal  mode  of  vibration.  This  process  is  called  resonance.  Here  it  is 
possible  to  get  large  structural  displacements  even  though  the  magnitude  of  the 
excitation  may  be  small.  If  the  magnitude  of  the  driving  force  is  large,  then  the 
structural  deformation  will  also  be  large.  It  is  also  important  whether  the  point  of 
application  of  the  driving  force  is  at  a  node  of  the  mode  shape  or  at  an  antinode.  If 
applied  at  a  node  the  response  is  less  than  if  the  force  was  applied  at  an  emtinode. 
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Figure  6-8.  ResiMnse  plot  of  the  barrel  tip  in  the  horizontal  direction 
(z-plane). 

Therefore,  finequency  response  functions  and  modal  analysis  are  crucial  in 
diagnosing  the  barrel  tip  displacement  of  the  PHALANX  gim. 

One  of  the  best  tools  for  understanding  which  normal  modes  are  involved  in  a 
forced  response  is  called  the  frequency  response  function  (FRF).  This  function  is 
obtained  by  applying  an  input  force  at  one  point  in  the  structure  and  measxiring 
the  response  such  as  displacement  at  another  point.  The  FRF  is  then  the  ratio 
the  response  divided  by  the  input.  In  this  study,  three  FRFs  were  obtained,  one 
for  response  of  the  barrel  tip  motion  along  each  coordinate  axis.  Mathematically, 
they  appear  as 

^3748^0})  Y 3748(Cfl)  (6  21) 

F,  (w)  ’  F,  (w)  ’  F,  (o)  ^  ' 
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In  these  three  functions,  the  numerator  stands  for  the  displacement  in  the 
respective  direction  at  node  number  3748  which  is  at  the  tip  of  the  bring  barrel. 
The  denominator  represents  the  driving  force  function  along  the  x  (fore-afl)  axis 
applied  at  node  number  1165  at  the  base  of  the  firing  barrel.  To  determine  the 
FRF,  the  firequency  of  the  force  is  swept  from  zero  to  a  high  value.  The  FRF 
appears  as  a  series  of  spikes  which  are  resonant  with  the  normal  modes  of 
vibration.  The  height  of  the  spikes  depends  on  the  magnitude  of  excitation  and  the 
width  of  the  peaks  depends  on  the  damping  associated  with  that  normal  mode.  A 
wider  spike  has  a  larger  damping  value.  This  FRF  is  complex  since  both 
amplitude  and  phase  of  oscillation  are  involved.  [Ref.  11] 

Figure  6-9  is  a  FRF  of  the  gun  response  along  the  x  direction  at  the  barrel 
tip.  As  expected  finom  the  transient  plot  in  Figure  6-6,  there  is  a  multimode 
response  involving  the  first  flexible  mode  at  14.4  Hz  and  second  at  15.5  Hz  which 
are  being  excited.  The  second  mode  is  much  more  excited  than  the  first  based  on 
the  height  of  the  two  peaks.  This  is  expected  since  the  second  mode  is  the  recoil 
mode,  while  the  first  is  a  rigid  body  pitching  mode.  It  appears  also  that  no  other 
modes  are  being  excited.  The  next  three  or  four  modes  involve  the  deformation  of 
the  gun  body,  and  do  not  appear  to  be  excited.  This  may  be  because  the  point  of 
appbcation  of  the  driving  force  is  at  the  rotor  which  is  suspended  by  the  recoil 
bearings. 

The  next  FRF  is  shown  in  Figure  6-10  describing  the  barrel  tip  response  in 
the  y  direction.  It  appears  to  verify  the  transient  response  in  that  this  is  a  single 
mode  excitation.  The  mode  that  appears  to  be  excited  is  the  recoil  mode  mixed 
with  some  pitching  motion.  This  is  identified  as  the  second  mode  due  to  the  fact 
that  its  normal  mode  motion  is  in  the  vertical  plane,  and  that  the  spike  is  fairly 
wide  suggesting  that  it  uses  the  higher  damping  value  of  0.33  rather  than  the 
others  at  0.2. 
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Figure  6-9.  Frequency  response  function  of  the  barrel  tip  in  the  x  direction. 

The  last  FRF  shown  in  Figiire  6-11  is  that  of  the  barrel  tip  in  the  z  direction. 
As  mentioned  earlier,  the  response  plot  of  Figure  6-8  showed  that  the  phenomenon 
of  beating  was  occurring.  Studying  the  FRF,  it  appears  that  the  first  flexible  mode 
at  14.4  Hz  and  third  flexible  mode  at  16.0  Hz  are  being  excited.  These  two  modes 
have  similar  damping  values  allowing  them  to  interfere  and  beat.  Beating  occurs 
when  two  modes  are  closely  spaced  in  firequency  and  amplitude  as  seen  here.  The 
amplitudes  were  different  in  the  x  direction  case,  so  that  beating  could  be  seen. 
The  two  modes  involved  in  the  z  direction  FRF  are  more  closely  excited  to  the  same 
amplitude  than  in  the  x  direction  where  only  one  mode  was  predominantly 
excited. 
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Figure  6-10.  Frequency  Response  Function  of  the  barrel  tip  in  the  y 
direction. 

F.  OTHER  FACTORS  AFFECTING  DISPERSION 

In  addition  to  the  resonance  described  above,  other  physical  phenomenon 
such  as  gyroscopic  whirl  and  rotor  imbalance  may  contribute  to  dispersion. 
Gyroscopic  whirl  is  a  self-excited  vibration  that  occurs  when  there  is  a  tangential 
force  normal  to  the  axial  direction  of  a  rotating  shaft.  At  some  critical  speed,  the 
whirl  instability  will  overcome  any  stabilizing  forces,  such  as  those  provided  by 
bearings.  The  amphtude  of  the  whirl  motion  will  only  be  limited  by  nonlinear 
stabihzing  forces.  This  type  of  motion  is  common  in  machinery  that  has  an 
overhanging  shaft  such  as  that  shown  in  Figure  6-12.  In  this  figure,  the  two  rear 
supports  can  be  thought  of  as  the  rear  needle  bearing  and  the  forward  thrust 
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Figure  6-11.  Frequency  Response  Function  of  the  barrel  tip  in  the  z 
direction. 

bearing.  The  mass  of  the  disc  at  the  end  of  the  shall  represents  the  gun  barrels 
and  clamps.  The  PHALANX  gun  system  is  similar  to  that  of  an  overhung  shall. 
This  fact  combined  with  the  knowledge  that  the  thrust  bearing  is  not  preloaded  are 
two  key  requirements  for  the  onset  of  whirling.  The  principal  excitation  force  of 
the  firing  gun  induces  several  moment  arms  which  would  lead  to  forces  normal  to 
the  axis  of  rotation.  These  indications  coupled  with  the  data  on  PHALANX 
dispersion,  such  as  the  pattern  shown  in  Figure  6-13,  strongly  implies  that 
whirling  is  occurring.  The  speed  of  the  whirl  is  usuedly  close  to  that  of  the  rotating 
shall,  and  accounts  for  the  fact  that  each  barrel  appears  to  shoot  in  the  same 
direction  on  each  cycle.  Discussions  with  Stan  Wilson  of  Port  Hueneme  indicates 
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that  field  experience  supports  this  possibihty  [Ref.  23].  When  an  old  thrust 
bearing  was  replaced  with  a  new  one,  the  dispersion  pattern  improved 
substantially.  This  would  seem  to  indicate  that  the  stabilizing  forces  of  a  new 
thrust  bearing  can  limit  the  onset  of  whirl  [Ref.  201 


Figure  6-12.  Representation  of  a  hung  shaft  demonstrating  the  concept 
of  gyroscopic  whirl 

Rotor  imbalance  is  another  phenomenon  found  in  rotating  machinery  which 
could  lead  to  barrel  tip  displacement.  If  the  principal  axis  of  rotation  does  not 
coincide  with  a  fixed  line  in  space,  then  wobbling  may  occur.  This  rotational 
instability  can  be  caused  by  machining  tolerances  or  mass  asymmetries.  Once 
started  it  can  aggravate  other  sources  of  vibration  when  rotational  speeds  are  high, 
as  in  the  PHALANX  gun.  The  unbalance  is  significant  since  it  determines  the 
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Figure  6-13.  A  dispersion  pattern  of  the  PHALANX  gun.  This  particular 
pattern  leads  to  assuming  that  whirl  is  occurring. 

degree  to  which  any  bending  or  flexural  mode  of  vibration  is  excited.  At  a  critical 
speed,  corresponding  to  a  resonant  frequency  of  the  system,  the  rotor  modal 
response  becomes  a  maximum.  When  a  rotor  spins  at  a  speed  near  its  critical 
speed,  it  tends  to  adopt  the  shape  of  the  normal  mode  that  is  excited.  The  amount 
of  rotor  displacement  depends  on  damping  in  the  system  as  well  as  the  modal 
components  of  the  imbalance.  This  modal  component  of  imbalance  may  be  reduced 
by  adding  a  nximber  of  correctional  tuning  masses  that  would  shift  the  normal 
mode  frequency  away  from  a  critical  speed. 
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Vn.  SUGGESTIONS  FOR  FUTURE  WORK 

A.  REFINING  THE  GUN  MODEL 

The  current  PHALANX  gun  model  is  a  first  cut  at  attempting  to  understand 
the  gun’s  complex  dynamic  behavior.  It  was  decided  early  that  it  was  better  to  go 
through  the  complete  simulation  cycle  from  finite  element  modeling  to  calculating 
its  dynamic  response,  rather  than  getting  caught  up  in  trying  to  make  a  perfect 
model  without  ever  calculating  its  response.  In  this  way,  the  whole  process  was 
validated  and  demonstrated  to  be  a  useful  tool  in  iinderstanding  the  dispersion  as 
well  as  modifying  various  components  to  see  how  the  real  gun  could  be  improved. 
Now  that  the  modeling  process  is  understood  finom  start  to  finish,  and  valuable 
experience  obtained  along  the  way,  certain  recommendations  stand  out  as  the  next 
logical  steps  in  continuing  this  research. 

The  first  recommendation  is  to  refine  the  gun  model.  The  purpose  of  this 
refinement  is  primarily  to  reduce  the  number  of  degrees  of  freedom  so  that  the 
model  would  have  a  smoother  animation  picture  and  a  quicker  solve  time.  The 
model  should  be  simplified  while  retaining  pertinent  characteristics.  This  would 
serve  to  better  illustrate  the  dynsmiic  behavior  of  the  gun’s  mode  shapes  and  lead 
to  more  insight  as  to  how  the  dispersion  pattern  is  being  formed.  Currently,  the 
animation  picture  is  slow  and  lurches  from  one  frame  to  the  next.  The  second 
purpose  behind  further  refinements  was  to  choose  better  elements  to  capture  even 
finer  details  of  the  gxin’s  behavior.  Now  that  the  real  gun  is  better  understood, 
some  components  of  the  gun  should  have  elements  that  model  bending  and 
torsional  behavior  with  beams  and  thin  shells  whereas  others  should  have  bricks  to 
model  elasticity  without  any  preferred  direction.  Neither  of  these  elements  can 
perform  both  functions.  Next,  the  gun’s  boundary  conditions  are  now  better 
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understood;  some  barrels  in  the  muzzle  clamp  are  constrained  while  others  are 
unconstrained.  Afler  attending  the  barrel  optimization  conference  in  Burlington, 
Vermont,  May  26th  and  27th,  1993,  the  consensus  was  to  firmly  clamp  all  the 
barrels.  This  modification  could  easily  be  evaluated  with  the  current  gun  model. 
Finally,  some  of  the  currently  modeled  components,  such  as  the  recoil  adapters, 
could  be  more  accurately  represented. 

1.  Muzzle  and  Mid-Barrel  Clamps 

In  an  effort  to  reduce  the  degrees  of  fireedom  of  the  model,  the  muzzle 
clamp  could  be  further  simplified.  Currently  the  muzzle  clamp  shown  in  Figure  3-9 
has  brick  elements  representing  the  plates  and  thin  shells  representing  the  webbed 
meshing  connecting  the  two.  This  arrangement  is  good  for  representing  the 
inertial  properties  of  the  clamp,  but  the  bricks  themselves  do  not  add  much  to  the 
model  other  than  increasing  its  complexity.  Since  it  is  unlikely  that  clamp  shear 
plays  a  significant  role  in  the  gun’s  dynamic  behavior,  such  elements  that  model 
shearing  well  are  not  necessary.  Bricks,  beams  and  shells  would  be  important  if 
this  model  was  used  to  understand  the  gun’s  thermodynamic  behavior,  but  that  is 
c\irrently  not  the  case.  Instead  it  is  proposed  that  the  bricks  be  replaced  by  a 
lumped  mass  located  at  the  center  of  gravity  of  the  muzzle  clamp  plates.  This 
lumped  maiss  is  an  "element"  choice  in  I-DEAS  where  a  structure’s  mass  and 
inertia  can  be  defined  at  a  geometric  point.  The  lumped  mass  can  then  be 
connected  to  the  barrels  via  rigid  elements.  Figure  7-1  shows  this  proposed 
arrangement.  These  rigid  elements  are  massless  and  do  not  flex,  but  they  are 
effective  in  attaching  two  structures  together.  In  this  case,  the  barrels  are 
attached  to  the  lumped  mass.  This  will  serve  the  purpose  of  defining  the 
appropriate  boundary  conditions.  Some  barrels  could  be  connected  while  others 
are  not  in  order  to  model  the  current  gun  in  the  fleet.  Alternately  they  all  could  be 
easily  connected  together  to  test  the  proposed  gun  design.  What  is  lost  in  this 
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representation  is  the  bending  behavior  of  the  plates.  It  is  uncertain  if  this 
behavior  is  important  to  model  or  not,  but  experiments  could  be  performed  on  the 
real  gun  at  the  Naval  Postgraduate  School  to  resolve  this  question.  In  a  separate 
finite  element  model  of  the  muzzle  clamp,  the  first  flexible  normad  mode  of 
vibration  occurred  at  about  500  Hz  and  is  far  fi-om  the  excitation  force  frequency  of 
75  Hz.  It  is  unlikely  that  this  mode  is  being  excited.  If  this  did  not  turn  out  to  be 
the  case  and  the  mode  is  excited,  thin  shells  could  be  used  to  represent  the  muzzle 
clamp  plates  vice  lumped  masses  and  rigid  elements.  Few  degrees  of  freedom 
would  be  added  to  the  model,  since  they  would  still  be  "cheaper"  to  use  than  bricks 
and  have  the  benefit  of  bending  behavior.  The  mid-barrel  clamp  would  be  modeled 
like  the  muzzle  clamp  plate. 

What  is  important  about  the  muzzle  and  mid-barrel  clamps  ^s  how  they 
determine  the  boxmdary  conditions  on  the  barrels.  Currently  three  of  the  barrels 
have  their  rotations  taken  out  at  the  muzzle  clamp  while  the  other  three  are 
allowed  to  rotate.  At  the  mid-barrel  clamp  the  barrels  are  rigidly  attached.  It  has 
been  proposed  that  the  PHALANX  gun  would  be  improved  if  its  barrels  were  all 
rigidly  constrained  at  the  muzzle  clamp  like  that  of  the  30  mm  GOALKEEPER  gvm 
[Ref.  21].  This  particular  gun  has  a  smaller  dispersion  pattern  than  the 
PHALANX  gun.  Therefore,  it  is  not  so  important  that  the  clamps  are  geometrically 
correct,  but  that  they  are  d)nnamically  correct. 

2.  Gun  Body  Housing 

The  next  area  of  finite  element  refining  is  that  of  the  gun  body  housing. 
Its  current  representation  can  be  seen  in  Figure  3-€.  The  model  could  be  improved 
if  there  was  better  matching  of  the  degrees  of  freedom  between  the  thin  shells  of 
the  gvm  body  and  the  bricks  of  the  thrust  bearing.  The  thin  shells  have  six  degrees 
of  fireedom  while  the  bricks  have  three.  This  creates  a  possible  boundary  condition 
problem  and  appears  to  manifest  itself  in  the  extreme  distortion  of  the  gun  body  in 
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Figure  7-1.  Recommended  modification  for  representing  the  muzzle 
I  clamp.  This  figure  shows  a  lumped  mass  connected  to 

the  barrels  via  rigid  elements. 

front  of  the  thrust  bearing.  Instead  of  obtaining  the  continuovis  d3mamic  behavior 

j 

of  the  gun  body,  it  has  been  broken  up  into  two  distinct  parts  with  the  thrust 
I  bearing  serving  as  the  divider.  The  forward  part  of  the  gun  body  acts  like  a  short 

dnim,  anchored  at  the  thrust  bearing.  By  itself,  its  finite  element  representation  is 
too  coarse  and  leads  to  unrealistic  modeUng  behavior.  This  behavior  appears  in 
the  higher  order  modes  of  vibration.  If  the  forward  gun  body  is  coupled  to  the  rear 
gun  body,  as  it  should  be,  then  its  finite  element  representation  is  not  too  coarse 
and  no  further  refinements  are  needed. 
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How  is  this  continuity  achieved?  The  brick  elements  used  to  represent  the 
thrust  bearing  should  be  replaced  by  beam  elements.  In  I-DEAS,  a  beam  element’s 
cross-section  can  be  defined  in  as  much  detail  as  desired.  The  inner  and  outer 
race,  as  shown  in  Figure  3-15,  could  be  defined  as  a  beam  cross-section.  By 
replacing  the  bricks  with  beams,  the  degrees  of  fi'eedom  are  then  matched  with  the 
thin  shells,  since  beams  also  have  six  degrees  of  freedom.  Therefore,  the  rear  gtm 
body  and  forward  gun  body  are  now  coupled  and  continuous  dynamic  behavior 
would  be  achieved. 

In  making  this  transition  to  beam  elements,  the  springs  used  to  model 
the  ball  bearings  would  also  have  to  be  modified.  Now  rigid  elements  would  be 
used  to  extend  firom  the  centroid  of  the  inner  and  outer  races  to  a  common, 
coincident  point  between  them.  The  springs  wotdd  be  co-located,  one  for  radial 
stiffness,  axial  stiffness  and  torsional  stiffiiess.  This  refinement  would  also  serve 
the  purpose  of  removing  any  moment  arms  that  currently  exist  in  the  model  due  to 
the  finite  distance  the  springs  extend  between  the  inner  and  outer  race. 

More  attention  could  be  used  when  modeling  the  details  of  the  gun  body 
housing.  Currently,  it  is  represented  as  a  cylinder,  whereas  in  reality,  there  are 
many  intricate  details  that  might  be  important  to  model  more  accurately.  For 
example,  the  ammunition  chuting  and  a  hydraulic  motor  coupled  to  the  gun  drive 
could  be  better  studied  now  that  the  gun  is  on-sight. 

3.  Rotor 

The  rotor  is  a  cylindrical  component  represented  with  brick  elements  as 
shown  in  Figure  3-16.  The  dynamics  of  the  rotor  are  not  critical  to  understanding 
the  gun,  but  could  be  better  represented  by  a  lumped  mass  and  inertia  properties 
while  keeping  a  portion  of  the  brick  elements  for  use  in  attaching  to  the  thrust  and 
needle  bearings.  In  this  way,  the  degrees  of  freedom  of  the  model  are  reduced, 
while  preserving  the  load  path  from  the  excitation  force  to  ground.  The  bricks 
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would  be  important  if  there  was  some  concern  that  shear  forces  were  causing  it  to 
deform,  or  if  thermodynamic  properties  were  being  modeled.  Since  neither  of  these 
characteristics  are  being  modeled,  the  use  of  brick  elements  for  the  whole  rotor  is 
unnecessary.  The  portions  used  to  couple  the  barrel  assembly  to  the  bearings 
should  be  kept,  and  should  be  connected  via  rigid  elements  to  a  lump)ed  mass  at  its 
center  of  gravity. 

4.  Stub  Rotor 

The  stub  rotor,  constructed  much  like  the  muzzle  clamp  as  Figure  3-11 
shows,  with  webbing  between  two  plates,  could  be  further  improved  by  dividing  the 
single  row  of  bricks  for  each  plate  into  two  rows.  The  stub  rotor  is  a  critical 
component  in  the  dynamics  of  the  gun,  perhaps  more  so  than  any  other  component. 
Its  rigidity  is  important  in  determining  how  much  the  gun  barrels  will  displace  and 
has  been  experimentally  studied  by  Martin  Marietta  to  calculate  its  characteristics 
[Ref.  16].  The  component  should  be  modeled  using  more  elements  instead  of  fewer 
elements.  By  including  two  rows  of  elements,  its  bending  behavior  can  be  better 
represented.  One  aspect  of  the  stub  rotor  and  barrel  clamps  that  is  difficult  to 
model  is  the  nonlinear  effects  that  gaps  play  in  the  gun  dynamics.  The  current 
version  of  I-DEAS  does  not  allow  this  option,  but  a  new  release  will  allow  the 
nonlinear  gaps  to  be  modeled. 

5.  Recoil  Adapters 

The  recoil  adapters  have  been  the  most  challenging  components  to  model. 
The  current  model  defined  them  in  terms  of  "ears"  that  stood  out  on  either  side  of 
the  gun  body  that  were  connected  to  groimd  via  rigid  elements  and  recoil  springs 
as  shown  in  Figure  3-16.  The  "ears"  were  solely  used  to  represent  the  adapter’s 
mass  and  to  model  the  correct  geometric  position  of  the  attachments.  The  model 
could  be  simplified  if  the  "ears"  were  replaced  with  lumped  mass  elements.  In 
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addition,  the  rigid  elements  could  connect  the  gun  body  to  the  springs  directly 
while  maintaining  the  correct  offset  of  the  pivot  points.  These  changes  are  needed 
in  order  to  provide  better  coupling  of  the  recoil  springs  to  the  gun  body.  As 
discussed  earlier,  problems  occur  when  there  is  a  discontinuity  between  degrees  of 
freedom  of  elements.  In  this  case,  it  occurs  between  the  thin  shells  of  the  body  to 
the  bricks  of  the  "ears."  This  recommendation  is  better  than  trying  to  accurately 
model  the  recoil  adapters.  It  is  likely  that  little  would  be  gained  for  such  detailed 
modeling  and  would  significantly  add  to  the  model’s  complexity. 

B.  VALIDATION  OF  THE  MODEL  THROUGH  TESTING 

After  constructing  the  finite  element  model  of  the  gun  and  obtaining  its 
normal  modes,  the  next  logical  step  would  be  to  compare  both  the  mode  shapes  and 
their  associated  frequency  with  those  of  the  real  gun.  However,  it  was  only  until 
late  in  this  study  that  a  M61A1  PHALANX  gun  was  received  at  the  school.  After 
proper  mounting,  a  modal  analysis  of  the  gun  could  then  be  performed  and  the 
results  compared  to  the  model.  The  equipment  already  exists  on-sight  for  this  type 
of  testing.  In  addition  to  testing  the  whole  gtm  as  one  unit,  it  could  be  broken 
down  into  a  barrel  cluster  and  gun  body  which  would  then  be  tested  separately.  In 
this  way,  discrepancies  could  be  localized  to  determine  which  components  needed 
to  be  better  modeled. 

Another  comparison  that  can  be  made  using  the  real  gun  is  that  of  finding  the 
proper  damping  values.  Although  the  damping  value  of  the  recoil  adapters  was 
provided  by  Martin  Marietta  and  verified  through  analytical  means  in  Chapter  VI, 
it  would  be  useful  to  obtain  these  values  experimentally.  It  is  uncertain  whether 
the  damping  values  vary  much  from  one  adapter  to  another,  or  the  effect  wear 
plays  on  them.  More  importantly  though,  all  other  normal  modes  of  the  gun  were 
estimated  to  have  a  damping  value  of  two  percent.  The  damping  values  of  the  gun 
model  could  be  correlated  with  the  real  gun  during  its  modal  analysis. 
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As  mentioned  earlier,  it  is  also  uncertain  how  much  free  play  exists  in  the 
gun  and  the  amount  each  component  contributes  to  the  total.  This  information 
could  also  be  determined  experimentally.  Individual  components  of  the  gun  should 
be  weighed  in  order  to  provide  another  means  of  comparison  so  that  the  modeling 
could  be  improved. 

C.  ADDITIONAL  MODELING 

1.  Modifying  the  Existing  Firing  Force 

In  an  efifort  to  get  through  a  complete  modeling  cycle,  there  were  many 
variations  of  the  existing  firing  force  that  were  put  off  until  a  later  time.  One  of 
the  modifications  which  would  certainly  add  to  the  accxnracy  of  the  model  is  to 
ramp  up  the  firing  force.  The  real  gun  does  not  immediately  jump  to  a  firing  rate 
of  4500  shots  per  minute.  This  could  be  accomplished  by  spacing  out  the  firing 
impulses  for  the  first  ten  to  twenty  shots  according  to  data  available  firom  the 
Naval  Surface  Warfare  Center  in  Dahlgren  or  Port  Hueneme.  Further 
modifications  would  be  to  model  the  diOerent  firing  rates  themselves.  This  study 
only  looked  at  the  firing  rate  of  4500  shots  per  minute.  However,  the  existing 
PHALANX  gun  can  select  a  rate  of  fire  of  either  3000  or  4500  shots  per  minute. 
There  has  been  a  proposal  to  investigate  a  firing  rate  of  6000  shots  per  minute.  If 
the  mass  of  the  penetrator  were  to  change  or  if  the  penetrator’s  velocity  were 
modified,  the  firing  force  would  change  and  could  easily  be  mod ''led.  The  gun  is 
capable  of  longer  engagements  and  longer  bursts  of  shots.  Longer  firing  bursts 
shoiild  be  modeled.  An  additional  and  perhaps  important  firing  force  to  model 
would  be  that  of  an  occasional  misfire.  In  this  case,  the  g\in  has  a  chance  to 
translate  forward  after  the  misfire  and  then  be  slammed  back  as  the  next  round  is 
fired. 
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2.  Adding  Additional  Forces 

In  this  study,  only  the  firing  force  was  represented.  The  impulse  of 
loading  each  round  may  also  affect  the  gun.  In  discussions  with  Martin  Marietta, 
this  force  is  not  insignificant  since  the  round  has  a  tight  fit  in  the  barrel.  This 
loading  force  is  applied  as  the  barrel  passes  through  its  lowest  point  of  circular 
motion.  While  the  fi:equency  of  the  loading  force  must  be  the  same  as  the  firing 
rate,  the  phase  relation  may  allow  it  to  add  to  or  subtract  from  the  firing  force. 

The  effects  of  rotation  can  also  be  modeled  in  I-DEAS.  Here  the  damping 
matrix  is  modified  to  represent  the  effect  of  rotation.  Even  though  the  modeled 
gun  is  not  rotating,  it  can  be  "tricked”  into  producing  the  same  dynamic  results  as 
if  it  were  rotating  by  modifying  its  damping  matrix.  This  involves  adding  off- 
diagonal  terms  in  the  viscous  damping  matrix.  The  normal  modes  of  vibration  are 
then  resolved  as  before,  and  the  excitation  forces  applied.  This  technique  is  only 
good  for  one  firequency,  and  the  model  will  have  to  be  resolved  for  each  different 
firing  rate. 

Another  force  to  model  is  the  variation  in  rotation  of  the  gim  by  adding 
an  impulse  perpendicular  to  its  axial  length.  Even  though  this  force  may  be  small, 
its  effect  could  be  large  since  it  might  be  a  contributing  perturbation  to  significant 
dynamic  forces  such  as  whirl. 

3.  Additional  Modeling 

With  the  current  PHALANX  gun  model  partially  validated,  it  is  possible 
to  move  on  and  use  it  as  a  tool.  The  U.  S.  Navy  is  interested  in  the  possibility  of 
reducing  the  PHALANX  gun’s  dispersion  by  constructing  a  gun  barrel  restraint  like 
the  one  shown  in  Figure  7-2.  This  structure  wiU  be  clamped  at  its  base  to  the 
isogrid  of  the  gun  mouint.  The  other  end  wiU  use  a  radial  bearing  to  clamp  onto 
the  gun  barrels.  Its  purpose  is  to  restrict  the  most  extreme  motions  of  the  barrel 
cluster.  A  finite  element  analysis  has  not  been  used  to  date  to  study  the  motion  of 
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thin  structure  and  optimize  its  design.  The  blueprints  have  been  received,  and  the 
Naval  Postgraduate  School  is  poised  to  take  on  this  task  [Ref.  24].  The  computer 
model  will  be  able  to  study  the  mode  shapes  of  this  combined  system  and  check  if 
there  is  any  problem  resonant  frequencies.  By  studying  the  mode  shapes,  it  wUl  be 
possible  to  determine  the  best  location  to  clamp  the  bearings.  This  system  can  be 
driven  and  damped  to  calculate  the  new  barrel  tip  displacements.  The  muzzle 
restraint  structure  could  be  modified  in  various  ways  to  find  its  optimum  design 
and  placement. 


for  Strength 


Figure  7-2.  Proposed  assembly  of  the  PHALANX  muzzle  restraint  [Ref.  24 


Another  modeling  task  is  to  swap  out  the  cxirrent  set  of  finite  element 
barrels  in  exchange  for  the  proposed  optimized  barrels.  This  will  help  determine 
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whether  the  new  barrels  will  pitch  more  than  the  old  ones.  The  new  barrels  will  be 
much  heavier,  nearly  18  kg  (39  lb)  compared  to  the  current  8  kg  (18  lb)  barrels. 
They  will  be  much  longer,  roughly  200  cm  (79  in)  compared  to  the  cxirrent  152  cm 
(59  in)  barrel  length.  The  effect  of  this  extra  weight  and  length  on  the  gun  mount 
and  its  dispersion  has  not  been  modeled  to  date.  The  NPS  model  would  be  an 
excellent  tool  for  making  this  analysis.  The  firing  force  would  have  to  be  modified 
to  take  into  account  the  extra  recoil  force  associated  with  the  heavier  enhanced- 
lethality  projectile  designed  to  be  used  by  these  barrels.  It  has  a  greater  mass  of 
0.290  kg  (0.64  lb)  than  the  old  projectile  mass  of  0.245  kg  (0.54  kg),  and  will  be 
fired  at  a  greater  velocity  increasing  the  specific  impulse  on  the  gun  mount.  An 
additional  modeling  step  would  be  to  use  the  proposed  muzzle  restraint,  as 
described  above,  with  these  new  barrels  and  solve  for  its  dynamic  response. 
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VIII.  CONCLUSIONS 


Using  a  finite  element  model  to  understand  the  complex  dynamics  of  the 
PHALANX  gun  has  made  possible  new  areas  of  analysis  and  component  design. 
This  has  been  the  first  effort  to  use  the  sophisticated  I-DEAS  computer  code  to 
study  the  gun.  According  to  the  Naval  Surface  Warfare  Centers  in  Dahlgren  and 
Port  Hueneme,  it  has  also  been  the  most  "technologically  advanced"  modeling 
research  of  this  important  weapon  system  [Ref.  25].  Not  only  does  it  provide  a 
thorough  and  detailed  model  of  the  gun  where  various  driving  forces  can  be  applied 
and  the  barrel  tip  response  calcxilated,  but  it  is  a  working  tool  to  assist  the 
PHALANX  program.  In  particvilar,  each  individual  component  of  the  gun  can  be 
studied  in  detail  to  understand  how  much  it  affects  the  overall  dispersion.  This  is 
an  especially  important  tool  today  when  budget  hmitations  allow  only  modest 
improvements  to  be  made  to  the  gun.  Furthermore,  new  components  can  be  tested 
and  their  performance  predicted  before  any  actual  hardware  tests.  This  can  be 
applied  to  the  proposed  muzzle  restraint  and  set  of  optimized  barrels,  before  any 
prototypes  are  ever  constructed.  The  cost  saving  and,  not  to  mention,  time  saving 
aspects  of  this  method  provides  an  invaluable  resource  to  the  program.  This  t}^ 
of  research  may  not  have  been  carried  out  before  by  private  industry  or  the  U.  S. 
Navy  because  the  computer  code  is  expensive  and  has  only  been  made  easily 
accessible  in  the  past  few  years.  The  Naval  Postgraduate  School  receives  an 
educational  discount  to  obtain  an  on-sight  hcense  for  the  code.  In  addition,  the 
high  concentration  of  professors  with  doctorates  and  U.  S.  Navy  officer/students 
makes  fur  a  unique  and  productive  environment.  The  school  has  been  proposed  as 
the  center  of  all  future  modeling  research  of  the  PHALANX  gun  [Ref  26].  In 
summary,  the  U.  S.  Navy  is  now  using  state-of-the-art  technology  to  understand 
and  improve  old  technology  in  a  cost  saving  way. 
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The  primary  goal  of  thia  study  has  been  to  construct  an  accurate  finite 
element  model  of  the  PHALANX  gun  and  demonstrate  the  potential  of  such  a 
sophisticated  computer  code.  This  goal  has  been  achieved.  Parts  of  the  model 
have  been  validated  with  the  real  gun,  and  the  results  of  this  study  have  been 
presented  to  the  several  organizations  involved  in  improving  the  PHALANX  gun 
[Ref  21].  The  model  is  poised  to  take  on  most  of  the  new  design  changes  that  the 
group  proposes.  The  model  is  now  considered  an  integral  part  of  the  decision 
making  process.  With  an  actual  g\in  on-sight,  the  model  can  be  completely  verified. 

Although  the  primary  purpose  of  constructing  this  finite  element  gun  model 
has  been  achieved,  some  preliminary  analysis  was  conducted  along  the  way.  One 
of  the  first  areas  of  concern  to  be  noticed  was  the  gun’s  normal  mode  of  vib'-ation 
around  14  Hz.  This  first  flexible  mode  with  associated  pitching  motion  would 
contribute  to  dispersion,  if  it  is  sufficiently  excited.  The  gun  rotates  at  a  frequency 
of  12.5  Hz  making  this  rotational  force  a  contender  for  exciting  this  mode. 
Frequency  response  functions  taken  of  the  barrel  tip  indicates  that  this  mode  is 
being  excited  by  the  principal  driving  force,  of  the  gun  firing.  This  discovery  was 
somewhat  unexpected,  but  upon  analysis  of  the  moment  arm  created  from  the 
driving  force  and  the  gun  moimting  positions,  it  is  reasonable  that  it  is  excited.  As 
the  g\ui  fires,  the  moment  arm  causes  the  gun  to  rock  back  and  pitch  upward 
coupling  the  firing  recoil  into  other  planes  of  motion.  Fourier  transforms  performed 
of  actual  dispersion  tests  appear  to  correlate  this  14  Hz  normal  mode  [Ref  22). 

As  far  as  design  analysis  of  the  g\in,  the  most  critical  component  contributing 
to  dispersion  was  revealed  by  this  study  as  the  double  angular  contact  bearing,  or 
thrust  bearing.  When  various  bearing  stiffiaess  values  were  tested,  the  effect  on 
barrel  tip  displacement  was  significant.  Confirmation  of  this  fact  came  from  the 
field  where  the  dispersion  pattern  improved  tremendously  as  a  worn  thrust  bearing 
was  replaced  with  a  new  one  [Ref.  23].  No  data  is  available  on  the  rate  of  wear  of 
the  thrust  bearing.  In  fact,  the  bearing  itself  is  not  inspected  while  on  ships,  but 
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only  when  the  gun  is  returned  for  depot  level  maintenance.  The  faci  that  the 
bearing  is  not  pre-loaded  is  a  leading  cause  of  the  physical  ejQTect  known  as  whirl 
that  occurs  on  an  overhung  shall  as  in  the  PHALANX  gun.  Dispersion  data  also 
suggests  the  possibihty  that  whirl  is  occurring.  The  addition  of  the  new,  heavier 
and  longer  barrels  may  significantly  aggravate  whirl  and  counteract  any  dispersion 
improvements  achieved  by  the  new  barrels,  especially  if  the  bearing  is  not 
redesigned  to  support  these  new  barreb.  Unfortunately,  this  appears  to  be  the 
case.  Not  only  would  the  added  torque  of  the  new  barrels  lead  to  a  greater  amount 
of  whirl,  but  it  would  probably  cause  a  greater  amount  of  wear  on  the  bearing 
which  would  further  worsen  the  problem. 

One  of  the  first  validation  tests  of  the  gun  model  was  to  drive  it  with  the 
principal  driving  force  of  the  rounds  firing.  In  the  response  plots  jhown  in  Chapter 
VI,  the  firing  gun  barrel  acted  as  predicted.  That  is,  the  gun  set  back  the  observed 
1.27  cm  and  oscillated  about  this  new  equilibrium  point.  Barrel  tip  displacement 
in  the  pitch  and  yaw  directions  corresponds  well  with  the  observed  2  milliradian 
dispersion.  Elliptical  dispersion  patterns  were  also  predicted  for  the  first  time. 

In  this  first  finite  element  study  of  the  PHALANX  gun,  the  potential  for  its 
future  use  has  been  documented  in  this  effort.  It  is  a  cost  effective  tool  that  has 
been  designed  to  help  the  U.  S.  Navy  not  only  understand  the  current  gun’s 
complex  dynamics,  but  those  of  any  future  g\in  as  well.  As  a  ship’s  last  layer  of 
defense,  the  PHALANX  gun  provides  a  valuable  function  for  the  ship  and  her  crew. 
This  is  especially  true  in  the  near  future  where  the  U.  S.  Navy  has  extended  its 
operating  envelope  into  new  environments,  where  the  l^t  layer  of  defense  may  be 
the  first  layer. 
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APPENDIX  A.  GUN  BARREL  NORMAL  MODES 


The  following  graphs  are  experimental  data  of  the  normal  modes  of  a  single 
PHALANX  gun  barrel.  The  barrel  was  suspended  by  two  low  stiffness  cords  and 
tapped  with  an  instrumented  hammer.  The  resulting  vibrational  information  was 
recorded  with  a  transducer  mounted  on  the  opposite  side  of  the  hammer  impact 
point.  The  normal  mode  frequencies  were  determined  through  3000  Hz  and  appear 
as  spikes  in  the  following  graphs. 
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APPENDIX  B.  GUN  MODEL  MODE  SHAPES 


The  following  pictures  are  snapshots  of  the  mode  shapes  used  in  the  dynamic 
analysis  of  the  gun.  There  are  20  flexible  modes  listed  in  ascending  order  by 
frequency. 
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